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Abstract. We study the one parameter family of Fredholm determinants det(J — 
Tocsin), 7 € R of an integrable Fredholm operator -ft" cs in acting on the interval (— s,s) 
whose kernel is a cubic generalization of the sine kernel which appears in random 
matrix theory. This Fredholm determinant appears in the description of the Fermi 
distribution of semiclassical non-equilibrium Fermi states in condensed matter physics 
^— ^ . as well as in random matrix theory. Using the Riemann-Hilbert method, we calculate 

' the large s-asymptotics of det (7 — 7-Kcsin) for all values of the real parameter 7. 
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1. Introduction 

2 

Let us consider the real vector space A4(n) ~ M. n of n x n Hermitian matrices 
M = (Mij) = M equipped with the probability distribution 



SI P^ N \M)dM = ce- NtrV ^dM, c J e- NtrV ^dM = l. (1.1) 

M(n) 

Here dM denotes the Haar measure on A4(n), N is a fixed integer and the potential 
V : R — > R is assumed to be real analytic satisfying the growth condition 

V(x) 

——2-— ^00 as|x|^oo. (1.2) 
[ — ■ m(x z + 1) 

00 ' 

The model (jVl(n),P( n ' N 'dM} is commonly referred to as unitary matrix model and it 
is well known (cf . \\.2\ [32] ) that basic statistical quantities such as 

O ■ E n ^(s) = Prob(Af G M.{n) has no eigenvalues in the interval (— s,s), s > 0) 

can be expressed as Fredholm determinants, indeed 

n . S S 

i"H ' ( 1)"^ f f 

; E n , N (s) = ^2 — ~ " det i. K ^N{xk, x^i ^dx! ■ ■ ■ dxj 

S : j=0 J ' s s 

= det(I-K n>N ) 
where K n ^ is the finite rank operator with kernel 

n— 1 „ 

K n , N (x,y) = e"T n*) e -% V(v) ^ Pl (x) Pl (y), / Pt (x) Pj (x) e - NV ^dx = 5 {j (1.3) 

i=0 £ 
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acting on L 2 ((— s, s), dX) . 

The given assumption (jl.2p enables one to study scaling limits n, N — > oo (compare 
[12], [16]). First the mean eigenvalue density ^K nt N(x,x) has a limit, that is for all 
bounded continuous functions g : R — > R 

lim - I g(x)K n N (x,x)dx = [ g{x)p v (x)dx = [ g(x)dp v (x) (1.4) 
n,N^Hx> n J J J 

and the support Xy of the equilibrium measure py consists of a finite union of intervals. 
Here the density pv(x) is determined by the potential V(x). At the same time, the local 
statistics of eigenvalues in the large n, N limit satisfies the so-called universality property, 
i.e. it is determined only by the local characteristics of the eigenvalue density py (see 
[371 El [III E]) : F° r instance let us choose a regular point x* G £y, i.e. py(x*) > 0. 
Then the bulk universality states that 

lim -A^Kn,n U + TTY' X * + ~ T^l = K ^ A = (1-5) 

rwoonpv^ar) \ npy{x*) npy[x*) J ir(A — p) 

uniformly on compact subsets of R, or in other words for a regular point x* 
lim Probl M G Ai(n) has no eigenvalues G (x* ■; — r,x* + 



\ V npy{x*) 1 npy(x*)- 

= det(I - K sin ) (1.6) 

where i^in is the trace class operator on L 2 ((— s, s); dA) with kernel K s i n (\, p) given 
in (|1.5p . The latter equation shows that in double scaling limits the basic statistical 
properties of hermitian random matrices are still expressible in terms of Fredholm deter- 
minants. On the other hand, the Fredholm determinant in the right hand side of (|1.6p 
admits the following asymptotic expansion [19] 



lndet(/-K sin ) = -^-iln(vr S ) + ^ln2 + 3C / (-l) + 0(^ 1 ), s -> oo, (1.7) 

where Q{z) is the Riemann zeta-function. This formula yields one of the most important 
results in random matrix theory, i.e. an explicit evaluation of the large gap probability, 
and its rigorous derivation was obtained in a series of papers [T4" l I2TJ 1 13"T]. 

The focus of the current paper lies on the computation of similar expansions for a 
certain family of generalized sine kernels. Let if C sm denote the trace class operator on 
L 2 ((— s, s); dX) with kernel 



sin 



'Hx 3 -p*) + x(x-p)) 



K csin (X,p) = ^ /x / — > xGM ( L8 ) 

7T(A - p) 

and det{I — K cs \ n ) the corresponding Fredholm determinant. This Fredholm determinant 
appeared in [9] in the asymptotical analysis of the Fredholm determinant 

det (I - K-pu) 

which corresponds to the first bulk critical universality class in the unitary matrix model: 
Let us assume py{x) vanishes quadratically at an interior point x* G Sy. Opposed to 
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(|1.5p . the scaling limit is more complicated OUT]. Let py(x*) = p' v (x*) = 0,p'y(x*) > 
and n, N — > oo such that 

lim n 2 / 3 (--l) =C 

n.,iV->oo \N J 

exists with C G R. Then the critical bulk universality guarantees existence of positive 
constants c and c\ such that 

lim — \j^K n \f ( x* H ^ 7 rr,3;*H ^-r-)=Kp U (X,p;x) (1.9) 

n,Af-i-oo cn 1 /^ \ cn L ' 6 cn 1 ' 6 J 

uniformly on compact subsets of R where the variable x is the scaling parameter defined 
by the relation 

lim re 2 / 3 ( 1 ) = xc\. 

n,N^oo \N ) 

Here the limiting kernel Kpu(X,p;x) is constructed out of the ^-function associated 
with a special solution of the second Painleve equation, the Hastings-McLeod solution 
u = u(x), see [9] for a precise definition of the Painleve II kernel. The main result in [9] 
is the following analogue of the Dyson formula (jl.T|> for the Painleve II determinant. As 

s — > oo 

oo 

2 13/" 
lndet(I — -ftTpn) = — -s 6 — s A x — -(sx) 2 — — Ins + (y — x)u 2 (y)dy 

X 

-^ln2 + 3C'(-l) + 0(^ 1 ), (1-10) 

and the error term in (jl.lUp is uniform on any compact subset of the set 

{i£l: -oo<x<oo}. (1-11) 

The proof of the latter expansion in [9] is based on a Riemann-Hilbert approach and an 
approximation argument which allows to derive the constant - with respect to s - term 

c = cq(x) = - lnF TW (x) - I In 2 + 3C'(-1) 

fa 

where Ftw( x ) is the celebrated Tracy- Widom distribution function, 

F TW {x) = e -ir(v-^ 2 (y) d y. 

In more detail, the x -dependence of the constant term follows from the analysis of the 
relevant differential identities. These identities though can not produce the numerical 
part of the constant. This part follows from the fact that in the large positive x-limit 

K Pn (X,p;x) = K csin (X, p; x) [ 1 + 0[x~ l/i e~i x3/ ) ) , x -> +oo, \,p£(-s,s). 



The last estimate, in conjunction with the Riemann-Hilbert analysis for the family K c 
with kernel 

- , s sin(|t(A 3 -/i 3 )+x(A-^)) r . 

K^ D (X,fJL;s) = ^ f j 1 ^, t G 0,1 

7r(A - p) 
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leads to the constant term in (jl.lOp . As a by-product of this analysis, [9] states, besides 
(jl.lOp . the asymptotic relation, 

lndet(/ - K csin ) = -\s & - s 4 x - \{sx) 2 - \ Ins - \ In 2 + 3C'(-1) + O(s^), (1.12) 
3 2 4 b 

valid as s — > oo and the error term in (|1.12p is uniform on any compact subset of the set 

In order to describe other spectral properties of large Hermitian matrices, one needs 
to study the Fredholm determinant 

det (I - 7 K Pll ) (1.13) 

for the values of 7 which are different from 7 = 1. Such one-parameter families of de- 
terminants already appear in connection with the sine- kernel determinant, for instance 
in the famous Montgomery-Odlyzko conjecture [35[ [36] concerning the zeros of the Rie- 
mann zeta-function, in the description of the emptiness formation probability and other 
correlation functions in one-dimensional impenetrable Bose gas [24[ [26] [27] as well as in 
a number of other important mathematical and theoretical physics applications. 

The analytical challenge of the determinants (|1.13p is once again the large s asymp- 
totics. In the case of the sine-kernel determinats, the result is well known (see |40[ [33j 
|34"1 [TJ [6] and [2] for more on the history of the question) 

1. As s — > 00, 

In det (I - jK sin ) = -Aktts + 2k 2 hi7rs + Xsin + 0(s~ 1 ) 
uniformly on and compact subset of the set {7 E R : —00 < 7 < 1}, where 

K = «(t) = ~^: ln ( 1 

The constant Xsin = Xsin (7) is given by the equation, 

P d 

Xsin = 2K 2 + 4K 2 ln2-4 / «(*)— argr(i«(t))cft, (1.14) 

Jo dt 

which was obtained by A. Budylin and V. Buslaev in [B] as a corollary of their 
main result in [6] - the asymptotics of the resolvent of the kernel -yK s [ n (X, /x). 
Formula (|1.14p also follows from the general theorem of E. Basor and H. Widom 
concerning the determinants of Toeplitz integral operators with piecewise con- 
tinuous symbols [2] . 

2. For 7 chosen from any compact subset of the set {7 G R : 1 < 7 < 00}, the 
Fredholm determinant det (7 — 7i^ s ; n ) has infinitely many zeros {s n } which ac- 
cumulate at infinity (see [3D ] [33" 1 13"4]). 

The main results of the present paper are the following analogues for the cubic sine- 
kernel (jl.8p . which together with f[l. 12[) state the large s behavior of det (J — jK cs i n ) for 
all values of the parameter 7. 
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Theorem 1. Let K cs \ n denote the trace class operator on L 2 ((— s, s); dXj with kernel as 
in OI). As s oo 



oo 

lndet(I-7ir csin ) = -k(^-s 3 + Axs^ +6k 2 \us- j (y-x)u 2 (y)dy + x + 0(s' 1 ) (1.15) 

X 

uniformly on any compact subset of the set 

{(7,x)eR 2 : -oo < 7 < 1, -oo < x < oo}, (1-16) 

where 

k = -w/(7) = -— ln(l - 7), 
P d 

X = 2k 2 + 8k 2 In 2- 4 / fs(t) — argr(ws(t))dt, 

with the Euler gamma-function T{z) and u = u(x,~f) denotes the real-valued Ablowitz- 
Segur solution of the second Painleve equation u xx = xu + 2u 3 corresponding to the 
monodromy surface 

M = {(si, . . . , s 6 )| si = -fy, s 2 = 0, s 3 = si, s n+3 = -s n }. 

The Ablowitz-Segur solution u = u(x,"f) [38] used in the statement of Theorem [1] is 
given as the unique solution to the boundary value problem 

u xx = xu + 2u 3 , u(x) ~ 7Ai(sc), x — > +00, 7 / 1. (1.17) 

Such solutions are smooth in case 7 < 1, with exponentially fast decay as x — > +00 and 
oscillatory behavior as x — > —00, see e.g. [22]. On the other hand for 7 > 1, the solution 
to (|1.17|) has poles on the real axis, but is still pole-free for sufficiently large positive x, 
in fact [1] for (7, x) chosen from any compact subset of the set 

|(7,x) G M 2 : 1< 7 < 00, x > Q ln 7^) j ( L1 8) 

the solution u = u(x, 7) to (|1.17p is pole-free. This in turn leads to our second asymptotic 
result. 

Theorem 2. For (7,2;) chosen from any compact subset of the set (|1.18|) . the Fredholm 
determinant det (J — 7-ftT C sin) has infinitely many zeros {s n } with asymptotic distribution 

80 1 , . . , / o \ r(l — in) 7r 

-s„ + 2xs n + - ln(7 - 1) In (l6s„ + 4xs n ) - arg v - ~ - + nvr, n 00. (1.19) 

3 7T '1 (ZKj 2 

The proofs of Theorem [1] and [2j are based on a Riemann-Hilbert approach. This 
approach (compare \25\ 115] ) uses the integrable form of the Fredholm operator i^ cs in) 
allowing us to connect the resolvent kernel to the solution of a Riemann-Hilbert problem. 
The latter can be analysed rigorously via the Deift-Zhou nonlinear steepest descent 
method. 

We should mention that a large class of the generalized sine-kernel determinants has 
already been considered in |30j (see eq.(1.6) there). In the case 7 < 1 and after a proper 
re-scaling, the determinant det (/ — 7-K"csin) can be put in the form which is very close to 
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the one treated in [30] . However, an essential difference occurs: the fast phase function, 
the function p(\) in the notations of [30] (see eq.(1.7)), which appear as a result of the re- 
scaling, does not satisfy one of the key conditions of [30]; moreover, it becomes depended 
on the large parameter. This means that the results of |30] are not directly applicable 
to our case. In fact, if one formally applies the main asymptotic formula of [30] to our 
case, then the first two terms of our asymptotic equation (|1.15p are reproduced while 
the constant (in s) term is not. Most significantly, the integral term with the Painleve 
function does not show up. The result of our Theorem[2j i.e. the asymptotic distribution 
(|1.19p of the zeros of the cubic sine-kernel determinant can not be compared with [30] 
since the techniques in [30] do not cover the situation 7 > 1. 

We also want to point out that the analysis of the cubic sine-kernel determinant 
is also of interest without the random matrix theory background: the determinant 
det (I — ifcsin) appears in condensed matter physics [5], namely in the description of 
the Fermi distribution of semiclassical non-equilibrium Fermi states. In order to under- 
stand perturbations to a degenerate Fermi gas one studies the one parameter extension 
of determinants corresponding to the kernel (|1.8p . that is 

det (I - 7 K csin ) , 7 6l. 

Our formula (|1.15j) provides therefore an exact formula for the large perturbations. 

Finally we would like to emphasize the appearance in our asymptotic results of the 
general Ablowitz-Segur solution of the second Painleve equation and not just the usual 
Hastings-McLeod solution. This fact, in particular, poses the following intriguing ques- 
tion which Theorem [2] does not answer: What is the large s asymptotic behavior of 
det (I — 7-fr cs i n ) when 7 > 1 and x coincides with one of the poles of the corresponding 
Ablowitz-Segur second Painleve transcendent ? 

Let us finish this introduction with a brief outline for the rest of the paper. Section 
[2] starts with a short review of the Riemann-Hilbert approach for the asymptotics of 
integrable Fredholm operators. We then apply the general framework to the Fredholm 
determinant det (/ — 7-fT cs in) and formulate the associated "master" Riemann-Hilbert 
problem (RHP). We also state logarithmic s,x and 7 derivatives of the determinant 
det (I — 7-Kcsin) and outline a derivation of an integrable system whose tau- function is 
represented by det (/ — jK cs [ n ). In sections I3 1151 following the Deift-Zhou roadmap, we 
construct the asymptotic solution of the master Riemann-Hilbert problem. At this point 
we will have to distinguish between 7 < 1 and 7 > 1, in fact for 7 > 1 we are lead to a 
solitonic Riemann-Hilbert problem. Using an extra "undressing" step, we overcome the 
singularities in the Riemann-Hilbert problem and derive the large zero distribution as 
stated in ([1.19]) . The situation is similar to the one dealt with in [10] . The calculations 
of sections [TU1 and [T71 provide us with the asymptotics of In det (/ — 7if C sin) including the 
constant term for 7 < 1. The large zero distribution ([1.19]) will be derived in section [TH1 

The authors dedicate this paper to the memory of Vladimir Savelievich Buslaev, whose 
pioneering works on the asymptotic analysis of integrable systems lie in the foundation 
of their modern asymptotic theory. 
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2. RlEMANN-HlLBERT APPROACH - SETUP AND REVIEW 



The given integral kernel (jl.8p belongs to an algebra of integrable operators first 
introduced in [25]: Let T be an oriented contour in the complex plane C such as a 
Jordan curve. We are interested in operators of the form XI + K on L 2 (T), where K 
denotes an integral operator with kernel 

g(A, M ) = £ ^ 1 . /<(A)W , y>(A)^(A) = 0, M G Z>i (2.1) 
A " ^ i=i 

with functions which are smooth up to the boundary of V. Given two operators 

XI + K , XI + K of this type, the composition (XI + K)(Xl + K) is again of the same 
form, hence we have a ring. Moreover let K l denote the real adjoint of K, i.e. 

A — fi 

Our results are based on the following observations (see e.g. [25] ). First an algebraic 
Lemma, showing that the resolvent of I — K is again integrable. 



Lemma 1. Given an operator I — K on L 2 (T) in the previous algebra with kernel (|2.ip . 
Suppose the inverse (I — K) -1 exists, then I + R = (I — K)~ l lies again in the same 
algebra with 

JH(A,M) = ££l fi A)gi(M) , EWi(A) = 



(2.2) 



and £/ie functions Fi,Hi are given by 

Fi(X) = ((J-^)-Vi)(A), fli(A) = ((/-^-^^(A). (2.3) 

Secondly an analytical Lemma, which connects integrable operators to a Riemann- 
Hilbert problem. 

Lemma 2. Ze£ -fC be of integrable type such that (I — K)~ l exists and let Y = Y(z) 
denote the unique solution of the following M x M Riemann-Hilbert problem (RHP) 

• Y(z) is analytic for z £ C\r 

• On the contour V , the boundary values of the function Y(z) satisfy the jump 
relation 

Y + (z) = Y_(z){l-27rif(z)h t (z)), zer 

where f(z) = (fi(z), . . . Jm{z)) and similarly h(z) = (hi(z), h M (z)) t 

• At an endpoint of the contour V, Y(z) has no more than a logarithmic singularity 

• As z — > oo 

Y(z) = I + 0(z~ 1 ) 
Then Y(z) determines the resolvent kernel via 

F(z) = Y(z)f(z), H(z) = (Y^z^hiz) (2.4) 
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and conversely the solution of the above RHP is expressible in terms of the function F(z) 
using the Cauchy integral 

Y(z)=I- [ F{w)h\w)-^-. (2.5) 
J w - z 

r 

Let us use this general setup in the given situation (|1.8|) . We have 

where y/z is defined on C\(— oo,0] with its branch fixed by the condition \fz > if 
z > 0. Lemma [2] provides us with the following K-RHP 

• Y(X) is analytic for A G C\[— s, s] 

• Along the line segment [— s,s), oriented from left to right, the following jump 
holds 

1 _ 7 7e 2i(|A 3 +xA) X 
_ 7e -2i(|A 3 +xA) j + 7 



y + (A) = y_(A)( _ 2i(| i +xA) r \ , ), A €[-*,*] 



• At the endpoints ±s, V(A) has logarithmic singularities, i.e. 

Y(X) = 0(ln(A=F*)), A -> ±a (2.6) 

• We have F(A) — >■ I as A — >■ oo. 
We can f actor ize the jump matrix 



- 7 e -2i(|A 3 +xA) ! + I e 3 ^ _ 7 i + 7 



and employ a first transformation. 

3. First transformation of the RHP 
We make the following substitution in the original A"-RHP 

X(\) = Y(X)e i{ l x3+xX)as , \eC\[-s,s]. (3.1) 
This leads to a RHP for the function X(X), our "master" RHP: 

• A"(A) is analytic for A G C\[— s,s] 

• The following jump holds 

X + (X) = X_(A) f 1 __~ 7 AG[- V ] (3.2) 



From ()2.6p . we deduce the following refined endpoint behavior 
X(X) = X(X) 



2m\-l I V A + s 



(3.3) 



where -A^(A) is analytic at A = ±s and the branch of the logarithm is fixed by 

X-s 
A+s 



the condition — ir < arg 4^ < ^- 
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• At infinity, X(X) is normalized as follows 

X(X) = (l+ 0(X' 1 )y i{ l x3+xX)a3 , X -> oo. (3.4) 

The latter master problem will be solved asymptotically as s — > oo in the next sec- 
tions by approximating its global solution with local model functions. Before we start 
this analysis in detail, we first connect the solution of lthe X-RHP to the Fredholm 
determinants det (/ — 7if C sin)- 

4. Logarithmic derivatives - connection to X-RHP 

We will express logarithmic derivatives of the determinant det(7 — jK cs [ n ) in terms of 
the solution of the X-RHP. To this end recall the following classical identity, valid for 
any differentiable family of trace class operators [39j 



ds 

In our situation 

dK, 



3 ( \ 

— In det (/ - 7-FTcsin) = -trace ( (l - Tocsin) ^ (7-^csin) j • (4.1) 

(A, fi) = K csin (X, y) (<5(// - s) + 6(ji + s)) , 



where, by definition 



/* 



8(w T s)f{w)dw = f(±s), 



and therefore 



-trace ((/ - jK csin ) 1 — ( / yK csia ) ) = -R(s, s) - R(-s, -s) 



with R(X, fj.) denoting the kernel (see (|2.2|> ) of the resolvent R = (I — ~yK cs i n ) ^ifcsin- 
The latter derivative can be simplified using the given definitions 

/i(A) = -h 2 (X), f 2 (X) = hi(X) 

as well as the identity detY(A) = 1, a direct consequence of the unimodularity of the 
jump matrix G(X) and Liouville's theorem: 

„ n , F l {X)H 1 ^) + F 2 (X)H 2 {^) F 1 (A)F 2 ( j »)-F 2 (A)F 1 ( / x) 

R{X,n) = = . 

A — fj, X — ii 

Since R(X,/i) is continuous along the diagonal A = fx (see ()2.2p ) we obtain further 

R(s,s) = F{(s)F 2 (s) - F^F^s), R(-s,-s) = F{(-s)F 2 (-s) - F^-s^-s) (4.2) 
provided Fi is analytic at A = ±s. One way to see this is a follows. Use (|3.ip and ([2 

F(X) = X(X)e- i ^ x3+xX ^f(X) 
to derive from ()3.3p the following local identity 



F(X) = 1(A) 



j(|A 3 +xA)(T3 

7 fi 



fW 



J W^i (43) 
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valid in a neighborhood of A = ±s. But this proves analyticity of F(X) at the endpoints 
and as we shall see later on, (|4,3p will allow us to connect (|4.ip via (|4.2|) to the solution 
of the X-RHP. We summarize 

Proposition 1. The logarithmic s-derivative of the given Fredholm determinant can be 
expressed as 

^lndet(/- 7 ^csin) = -R(s,s)- R(-s,-s), (4.4) 

R{±s,±s) = F[(±s)F 2 (±s) - F^±s)F 1 (±s) 
where the connection to the X-RHP is established through 



lndet(7 - jK csin ) = -trace f(l - jK csia ) (jK csin ) ^ . 



Besides the logarithmic s-derivative we also differentiate with respect to x 
d_ 

dx 

In our situation the kernel depends explicitly on x, indeed 

|-( 7 K csin ( A ,M)) = i(/i(A)/n(A0 " / 2 (A)Mm)) 

and with ([23]) 

s 

-trace {{I - 7 X csin )~ 1 ^^csin) = -i f (F 1 (X)h 1 (X) - F 2 (X)h 2 (X))dX. 

—s 

On the other hand the Cauchy integral (|2.4[) implies 

s 

y(A) = / + ^ + 0(A~ 2 ), A^oo;mi = j F{w)h t {w)dw 

—s 

SO 

9 

— lndet (I - 7-fCcsin) = i{m\ - mj 1 ), mi = (m\ J ) 
and the connection to the X-RHP is established via (|3.4p . Again we summarize 



Proposition 2. 27ie logarithmic x-derivative of the given Fredholm determinant can be 
expressed as 

^ lndet(J - 7 ^csm) = i(Xf 2 - X?) (4.5) 

with 

X(A)~(/ + ^ + § + § + 0(A- 4 ))e^ 3 +^, A^oo; X 1 = (X»). 



A A 2 A 3 

The logarithmic s and x derivatives are sufficient to determine the large s-asymptotics 
of det(I — 7i^ cs ; n ) up to the constant term. In order to evaluate the constant explicitly 
for 7 < 1, we use the following approach (see [13] for a similar method used in the 
asymptotics of Toeplitz determinants). 
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s 

^- lndet(J - 7-PCcsin) = -trace - Tocsin) _1 -Kcsin J = ~ J R(K A)dA. 



Our goal is to express the latter integral over the resolvent kernel in terms of the solution 
of the underlying X-RHP. Recall (gS D , (g3D , the definition of the functions f(X),h(X) 
and unimodularity of Y(\) 



R(X,X) = F{(X)F 2 (X)-F^(X)F 1 (X) = ^{AX 2 + x) + ^-{Y{ 1 (X)Y 22 (X) 

-y n (A)y 2 ' 2 (A) + y/ 2 (A)y 21 (A) - y^(\)y 12 (x)) + (y{ 1 (x)y 21 (x) 
-y 11 (A)y 2 ' 1 (A))/ 1 2 (A) + (y/ 2 (A)y 22 (A) - yi 2 (A)y 2 ' 2 (A))/ 2 2 (A) 

where (') indicates differentiation with respect to A. Next (|3.ip and unimodularity allows 
us to rewrite the previous expression for R(X, A) solely in terms of the X-RHP 



fl(A,A) = ^,[{X' n (X)+X[ 2 (X)){X 21 (X) + X 22 (X)) 
-{X n (X) + X 12 (X)){X' 21 (X) + X 22 (X)) 



(4.6) 



Our next move will replace all terms involving derivatives with respect to A and this can 
be done by considering the differential equations associated with the X-RHP, see also 
section [5] below: 

All jump matrices in the X-RHP are unimodular and constant with respect to A, thus 
the well-defined logarithmic derivative X^X _1 (A) are rational functions. Indeed using 
(ED and (GOD 



dX 
dX 



, b\ , (d e A A B 



X = A(X,s,x)X (4.7) 



where 



—X(s) 



and with parameters b, c, d, e, f which can be expressed in terms of the entries of X\ and 
X 2 



b = 2X\ 2 , c = 2Xf , d = ix + UX\ 2 X 21 



(4.9) 



e = 8i(xl 2 Xf 2 - X 2 12 ) , / = Sifx^Xl 1 - X 2 21 ) . (4.10) 
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Substituting (|27fj) into P~6j) and recalling (gjHD we obtain with A = (Aij),B = (By) 



^(A,A) = 

27rz 



8iA 2 + 2d + 



Au — A22 B\\ — B22 



X-s 



X + s 



(X n (X)+X 12 (X)) 



(X 2 i(X) + X 22 (X)) 
+ ( - AiXb + e + 



A\2 B\2 \ tv />x , ^ ^\\2 



A — s X + s 



A — s X + s 



(x 21 (x) + x 2 2(x)y 

(X n (X) + X 12 (X)f 



(4.11) 



Next we 7-differentiate the X-RHP in (|3.2[) to obtain the following additive RHP for 
the function H(X) = ||(A)(X(A)) _1 

• H(X) is analytic for A G C\[— s, s] 

• Along the line segment [— s, s], oriented from left to right 



1 1 



(A)=fl-_(A) + X_(A)f_J J (X_(A))"\ AG 



-s, s 



H(X) has at most logarithmic singularities at the endpoints A = +s 



dX 



H(X) = —(X)(X(X))-' + X(X)^ 
• As A — > 00, we have H(X) — > 



1 f-1 1 



X-s 



J Jj (X(A)) In^, A^±s (4.12) 



If we let 



1 1 



K(X) = X(X)[ 1 1 \(X(X)) \ XeC\[-s,s], 



then K + (X) = K_(X), X G [— s, s] and K(X) is bounded as A — > +s. Hence K(X) is entire 
and we have a solution to the if -RHP 



1 f K_(w) 



w — X 



dw 



1 f K(w) 



2iri I w — X 



dw 



1 

2^1 



(Xn(«;)+A'i2(to)) 2 



-(X 11 (w)+X 12 (w))(X 21 (w)+X 22 (w)) 

-(X 21 (w)+X 22 (w)) 2 (X 11 (w)+X 12 (w))(X 21 (w)+X 22 (w)) I w - \' 



dw 



This solution enables us to rewrite J^ s R(X, X)dX in (|4.1ip . for instance 

A n (X 11 (A) + X 12 (A))(X 21 (A) + X 22 (A))dA= I w n H ll (w)dw, n G Z> 



with £ denoting a closed Jordan curve around the interval [— s, s] and where we used 



1 



A n = — 



2iri / w — X 



dw, X G [— s, s]. 
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Similarly 



X n (X 21 (X)+X 22 (X)) 2 dX 
X n {X u (X)+X 12 (X)) 2 dX 



w n H 2 \(w)dw, 
j w n H 12 (w)d 



and we obtain 



d_ 



1 

'27d 



8i J w 2 Hn(w)dw 



lndet(/- 7 ^csin) = -- / R(X,X)dX 
7 J 

— s 

+ J (2dHn(w) + eH 21 (w) + fH 12 (w))dw - 4i J w(bH 21 (w) - cH 12 (w))dw 
s s 

s 

- J {[A n - A 22 )K U (X) + A 12 K 21 (X) + A 2l K l2 (X) 

—s 
s 

+ J [(B n - B 22 )K U (X) + B 12 K 21 (X) + B 21 K 12 (X) 



dX 
X-s 

dX 



X + s 



(4.13) 



Since 



and 



27T? 2vT7 

K (A) = — A + 0(A - s), A ->■ s, if(A) = — 5 + OfA + s), A ->■ -s 

7 7 



(A u - A 22 )A U + 2,4 12 A 21 = = (B u - B 22 )B U + 2B 12 B 21 , 



we deduce that the last two integrals in (|4.13p are indeed well-defined. To evaluate them, 
let 

H(X) = H(X) - X(s) (Z\ jW^JLjn^, XeC\[-s,s] 



H(X) = H{X)-X(-s) 



2iri X + s 



From ()4.12|) we see that i?(A) is bounded as A — > s and H(X) is bounded as A — > —s, 
more precisely 

H(s) = ^(s){X(s))-\ H(-s) = ^(- s )(X(-s)y\ 



also 



H + (X)=H-(X)+K(X)-K(s), H + (X)=H-(X)+K(X)-K(-s), AG [-a, a], 
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hence 

^■IfaiA g(A)-- L/^-y-)*, (4.14) 
27ri j u; — A 2iri J w — X 

—s —s 

and we conclude 

s 

dX 



J ((A n - A 22 )K U (X) + A l2 K 2l {X) + A 21 K 12 (A) 



X-s 



7 A - s y_ s A - s 

+^2i j Kl2(X) x ~_ j 12(S) = 2vri [(A n - A 22 )H u (s) + A 12 ff 21 (s) + A 21 ff 12 (s)] . 

— s 

Similarly 

s 

y ((Bn - B 22 )K U (X) + B 12 K 21 (X) + S 21 K 12 (A))^- 

= 2™ [(flu - B 22 )H u (-s) + B 12 H 21 (-s) + fl 2 ifli 2 (-s)] . 
To evaluate the remaining integrals in (|4.13p . we recall 



l,„ t 1 



H(X) = T (X 1 ) 7 + - 5 ((X 2 ) 7 -(Jf 1 ) 7 X 



A v 1/7 ' A 2 

+ T3 ((^3) 7 + (*i) 7 (*i - * 2 ) - (^2)7^1) + 0{X^), X ->• 00 



A 3 

and apply residue theorem 

y H(w)dw = 2vri(Xi) 7 , y wH{w)dw = 2-ni((X 2 ) 1 - (Xi) 7 X x ) 
s s 

y w 2 H{w)dw = 27ri((X 3 ) 7 + (X 1 ) 7 (X 1 2 - X 2 ) - (X 2 ) 7 Xi). 
s 

We summarize the previous computations 

Proposition 3. The logarithmic ^-derivative of the given Fredholm determinant can be 
expressed as 

F) 11 
— lndet (/ - 7 ^csin) = -8i((X 3 ) 7 + (^i) 7 (^i 2 - X 2 ) - (X 2 ) 7 Xi) 

-2d((X 1 ) 7 ) U + 4ib((X 2 ) 7 - (X!) 7 X X ) 21 - 4ic((X 2 ) 7 - (X!) 7 Xi) 12 

-e((X!) 7 ) 21 - /((Xi) 7 ) 12 + ((An - A 22 )H u (s) + A 12 fl 21 (s) + A 2l H l2 (s)) 

(B u - B 22 )H u (-s) + B 12 H 21 (-s) + B 21 H 12 (-s)) (4.15) 
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with 



X(X) 



oo 



A ' A 2 A 3 

and the functions b,c,d,e, f, A, B are defined in (|4.9|) . (|4.10|) and 

5. Differential equations associated with the determinant det (I - 7-Ff C sm) 

Our considerations rely only on the underlying Riemann-Hilbert problem. Neverthe- 
less, before we continue the asymptotical analysis, we will take a short look into the 
differential equations associated with the X-RHP. 

We already used the differential equation 

dX 

— = A(X, s, x)X 

in the derivation of (|4.15p . Also, the logarithmic derivatives X S X~ 1 (X) and X X X~ 1 (X) 
are rational functions, indeed 



and furthermore 



dX 

ds 

dX 

dx 



A 



B 



A + s 



iXa? 



b 
-c 



X = Q(X,s,x)X 



X = U{X,s,x)X. 



Hence we arrive at the Lax-system for the function X 

dx 
dX 



• ^ = A(X,s,x)X 



dX 
ds 

dX 



@(X,s,x)X, 



■jz -U{X,s,x)X. 

Considering the compatibility conditions 

A s -e x = [&,A], A X -U X = [U,A], @ x — U s = [n, <d] (5.1) 

we are lead to a system of eighteen nonlinear ordinary differential equations. Since it is 
possible to express the previous derivatives of In det (I — 7i^ cs ; n ) solely in terms of the 
unknowns b, c, d, e, /, A and B. Since it is possible to express the previous derivatives of 
In det (/ — 7-KTcsin) solely in terms of the unknowns b, c, d, e, /, A and B one can then try 
to derive a differential equation for the Fredholm determinant using (|5,ip . We plan to 
address this issue in a future publication. 

6. Second transformation of the RHP - rescaling and opening of lenses 
Let us scale the variables in (|3.1|) as A = zs, so that 

T(z) = X(zs)e- sH ^\ z€C\[-l,l], = i (^z 3 + ^\ (6.1) 

solves the following RHP, which up to the rescaling A = zs, is identical to the initial 

y-RHP 

• T(z) is analytic for z G C\[— 1, 1] 
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Along the line segment [—1,1] oriented from left to right 



T + {z) = T^z)e sH ^( 1 - 1 , 7 

\ -7 1 + 7 

In a neighborhood of the endpoints z = ±1 



*€ [-1,1] 



T( z )e sH ^ = X{zs) 
At infinity 



I 



7 

2iri 



-1 1 
-1 1 



In 



z - 1 
z+l 



T(z) = I + 0(z~ 1 ), z^oo 



(6.2) 



We now move to a RHP formulated according to the sign-diagrm of Re which is 
depicted in Figure [TJ In this Figure we choose x from a compact subset of the real line, 
s > is sufficiently large and 



V3 v 

denote the two vertices of the depicted algebraic curves. With the matrix factorization 




Figure 1. Sign-diagram for the function Re In the left picture 

we indicate the location of z± as x > and in the right picture for a 
particular choice of x < 0. Along the solid lines Re $(z) = and the 
dashed lines resemble arg z = ±^ 



-7 7 A ( 1 0V A 7(1-7)" 

-7 1 + 7^ HI"?)- 1 V 1 7j V° 1 

= SlSdSu, (6.3) 
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valid as long as 7 7^ 1, we perform opening of lenses as follows. Let L ■ denote the upper 
(lower) lense, shown in Figure El which is bounded by the contours 7 - fe with 




732 741 



Figure 2. The second transformation - opening of lenses 



S(z)L(z) 



Define 

' T(z)S^\ z g £+u£+ 
S{z) = I T(z)S L , zeCzULj, 
T(z), otherwise, 
then S(z) solves the following RHP 

• S(z) is analytic for z £ C\([-l, 1] U C), with C = \Ji,j {iji U 7^) 

• The following jumps hold, with orientation fixed as in Figure [2] 

S+(z) = S^{z)e sH{z)rj:i G s e- sH{z)cT \ z G C\([-l, 1] U C) 
where the piecewise constant matrix Gs is given by 



G s 



q-l 

Su, 

a-7r 3 , 

Sl, 
0-1 



2 G 7n u 72i> 
z G 7^ U7 2 + 2 , 

ze [-1,1] 

« G 7ii U7I1, 

Z G 732 U 742 • 



• As 2 — > ±1, we have 
S(z)ir x {z)e 



X(zs) 



1 + 



At infinity, S(z) = 1 + 0(z~ l ), 



7 

2vri 
00. 



■1 1 
■1 1 



In 



z+1 



(6.4) 



■5) 



(6.6) 



ASYMPTOTICS OF A CUBIC SINE KERNEL DETERMINANT 18 

Let us analyse the behavior of the jump matrix Gs(z) in (j6.5[) along the infinite branches 
as s — > oo. To this end recall the sign-diagram of the function Rei)(z), depicted in Figure 
[TJ We have in the upper half-plane 

G 5 (z) = e s3 ^ 3 (j jV^*)"*, zG 7 J, j,k = 1,2 

with a constant a G C. Since we choose a; from a compact subset of the real line and 
s > is sufficiently large, Rei?(z) is always negatice on 7^,,j, A; = 1,2 outside small 
neighborhoods around the origin and the endpoints z = ±1. Hence for such z 

G s (z) — >I, s^oo (6.7) 

uniformly on any compact subset of the set (|1.16p and the stated convergence is in fact 
exponentially fast. As similar statement holds on the infinite branches in the lower 
half-plane. There 

G s (z) = e^^(l J)e-^)-, z€ ^ j, k = l,2 

again with some constant b G C. In this situation Re$(z) > outside a small neighbor- 
hood of the origin and therefore 

G s {z) — ► /, s -)■ oo (6.8) 

also uniformly on any compact subset of the set (|1.16p . From (|6.7p and (|6.8p we expect, 
and this will be justified rigorously, that as s — > oo, S(z) converges to a solution of the 
model RHP, in which we only have to deal with the diagonal jump matrix Sd on the 
line segment [—1,1]. Let us now consider this model problem. 

7. The model RHP 

Find the piecewise analytic 2x2 matrix valued function M(z) such that 

• M{z) is analytic for C\[-l, 1] 

• The following jump condition holds along [—1, 1] 

M+(z) =M-(z)S D , zg[-1,1] 

where 

s D = (i- 7 r 

• M(z) has a most logarithmic singularities at the endpoints z = ±1 

• As z — > oo we have M(z) = I + 0(z _1 ) , z — > oo 

Only assuming that 7 / 1 we can always solve this diagonal and thus quasi-scalar RHP 
at once (see [22]) 

1 



M(z) = exp 



1 f ln(l - 7) CT3 



2iri I w 



dw 



z + 1 
z-1 



-1/0-3 



(7.1) 



v=— ln(l-7), arg (1 - 7) G (-7T, vr] 

and is defined on C\[— 1, 1] with its branch fixed by the condition (fr^)" — > 1 as 

z — > 00. 
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Remark 1. We bring the reader's attention to the important fact, that in case 7 < 1, 
we have arg(l — 7) = and v is therefore purely imaginary. However if 7 > 1, then 
arg (1 — 7) = 7r and v equals 

1 , / x 1 1 

^ = ^— In ( 7 -l) + - = i/ +2> i/ ° e *' 

Later on we will see that this difference will have a very substantial impact on the whole 
steepest descent analysis. 

8. Construction of a parametrix at the origin z = 

In this section we construct a parametrix for the origin z = 0. This construction 
involves the Ablowitz-Segur solution u = u(x, 7) of the boundary value problem 

u xx = xu + 2u 3 , u(x) ~ 7Ai(x), x — > +00, 

where Ai(x) is the Airy-function. Viewing x,u and u x = dU jL as real parameters, 
consider the 2x2 system of linear ordinary differential equations, 

a /j- \ t At j- \ a -a2 , -a / ?A f—ix — 2iu 2 —2u x \ 
— = A( C , x)*, A( C> = -4< 2 a 3 +4( (_ u Q j + ( ix + ^ J • (8.1) 

The given system has precisely one irregular singular point of Poincare rank 3 at infinity, 
thus from the classical theory of ordinary differential equations in the complex plane, 
we obtain existence of seven canonical solutions ^ n (C) which are fixed uniquely by their 
asymptotics (for more detail see [22] ) 

* n (C)~ (/ + 0(r 1 ))e-^ 3+ <^, C^oo, (£!l n 

in the canonical sectors fi n given in Figure [3] 

fl n = {CeC| argCG (|(n-2),|n),n = l,...,7}. 




Figure 3. Canonical sectors of system (|8.ip with the dashed lines indi- 
cating where Re £ 3 = 
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Moreover the presence of an irregular singularity gives us a non-trivial Stokes phe- 
nomenon described by the Stokes matrices S n : 

S n = (^(C)) _1 *n+l(C)- (8-2) 

In the given situation (|8.ip with the Ablowitz-Segur solution u = u(x,j), these multi- 
pliers are 

Si = ij ; ^2 = jj , S4 = , S3 = Si, S5 = S2, S$ = §4. (8.3) 

We now choose Pn(C) to be the first canonical solution of ()8.ip corresponding to the 
specific choice (|8.3|> . i.e. 

P//(C) = *l(C), argCG (~ I 



and assemble a piecewise analytic model function as follows 

PiT(C) = < 



Pn(C), argCG(-f,f)U(f,^), 



P//(C)<Si, argCG(f,f), (8.4) 

77T ll7T\ 



{ P//(C)5 4 , argCG(f, i f 
One checks directly that Pjj(C) solves the model RHP shown in Figured] 




Figure 4. The model RHP near z = which can be solved explicitly 
using the real-valued Ablowitz-Segur solution of the second Painleve 
equation 

• p7P(C) is analytic for ( € C\{arg C = f , f , 7 f , ^} 

• On the infinite rays, the following jumps hold 



(P/f (C)) + 


= (PjT(O) 


_Si, 


arg 


;C = 


7T 

6 


(PiT(0) + 


= (PiTiO) 


St 1 , 


arg 


C = 


T 


(PiT(0) + 


= (PiTiO) 


_Si, 


arg 


c = 




(PiT(0) + 


= (PiTiO) 


_°4 ; 


arg 


c = 


llvr 
IT 
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In terms of the previous choice (|8.3p . we have the following uniform asymptotics, 
valid in a full neighborhood of infinity 

^(0~(/ + ^ + ^ + ^ + O(r 4 ))e-^ 3+ ^, C^oo (8.5) 

with 

1 f—iv u\ If u 2 — v 2 2i(u x + uv 

iiu\ — — 

and 



7711 " 2 1 « t" 7 ' m2 " 8 I "2*K + uv) u 2 - v 2 



1 / i(v 3 - 3vu 2 + 2(xv - uu x )) -3(u(u 2 + v 2 ) + 2(vu x + xu)) 
18 V-3(u(n 2 +^ 

where we introduced 



7713 48 1 -3(ii(u 2 + v 2 ) + 2(mi a! + xu)) -i(v 3 - 3vu 2 + 2(xv - uu x )) 



l 



V = (^a;) 2 — XU 2 — U 

Remark 2. In view of our discussion in section {J\ the Ablowitz-Segur solution u = 
u(x,7) might have poles on the real line. However this can only happen in case 7 > 1 and 
in this situation we restrict ourselves to values of (7, x) from (|1.18p . Thus in either of our 
cases, u = u(x,j) is smooth in x and therefore the model function Pj] H '(C) = Pn H {C'i x ) 
well defined. 



The model function (|8.4|) will now be used to construct the parametrix to the solution 
of the original S'-RHP in a neighborhood of z = 0. First set 



PiT(C) 



e *iu<T 3 pKH (Q e -™°a ^ ImC>0; 
^ e TTivo z pRH (£) e irii/os a ImC<0; 

where v was introduced in (|7.ip . This leads to a model RHP with jumps on the positive 
oriented real line 

(PiT(0) + = (pl t i H (())Ai-ir a:i 

as well as on the infinite rays arg C = 5 ' 



6 ' 6 ' 



(^(C)) + 


= (p?i h (0) 


(^(C)) + 


= (pFK)) 


(PiT(0) + 


= (PiTiO) 


(PiT(0) + 


= (p?i h (0) 



1 0\ _ 7T 

-^(l-T)- 1 V' argC "6 

1 (A . 5vr 

arg C 



1 i7(l-7) _1 \ A 77r 

arg C = — 



1 J ' b s 6 

1 - !7(1 r 7,_1 ). -c=!£ 



and with behavior at infinity 



77 lU V + C C 2 C 3 1 V I e 2TO , ImC<0; 
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where 



1 / -iv ue 2lTiu \ _ If u 2 - v 2 2i(u x + uv)e 27Tiu 

in — - 

and 



mi ~ 2 \ue- 27Tiu iv h 1712 ~ 8 \ -2i(u x + uv)e~ 2 ™ u 2 - w 2 



J_ / i(v 3 - 3vu 2 + 2(aiv - uu x )) -3(u(u 2 + t; 2 ) + 2(vu x + xii))e 2?ri,/ 

m3 ~ 48 \-3(u(u 2 + v 2 ) + 2{vu x + xujje - ** -i(^ 3 - 3vu 2 + 2{xv - uu x )) 

Secondly define 

CO) = sz, \z\ < r (8.6) 

which yields a locally conformal change of variables and which allows us to define the 
parametrix U (z) near z = by the formula 

^(«) = ffiS (z)e- < 5^i^(C(«))e i <s c W + * c W^e i T^<7 1 , \z\ < r, a x = K J J (8.7) 

with C(^) as in fj8.6j) and the matrix multiplier 

By construction, in particular since A)( z ) is analytic in a neighborhood of z = 0, the 
parametrix U(z) has jumps along the curves depicted in Figure [5J and we can always 
locally match the latter curves with the jump curves of the original RHP. Also these 
jumps are described by the same Stokes matrices as in the original S'-RHP, since the 
previously stated jumps of Pj\ H (C) w iU be conjug ated with e^ " 3 aie~ sH ^ which pre- 
cisely matches model jumps with original ones. Thus the ratio of S{z) with U(z) is 
locally analytic, i.e. 

S(z) = N (z)U(z), \z\ < r < X - 





Figure 5. Jump graph of the parametrix U(z) 



Let us explain the role of the left multiplier Bq(z) in the definition (|8.7p . Observe 
that 

I, Im C > 0; 



ffiB («)ffi = M(z) 



e 2™<7 3) l m ^<0. 
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This relation together with the asymptotic condition (|8 . 5[) implies that, 

U(z) = a l B Q {z)e-^(l+ 1 ^ + ^ + |* + 0(C 4 ))e^ B^ 1 (z)a 1 M(z) 



t . i j3 ( \—lf v ue~ 2 ™\ 

I + -B (z) l^_ ue2mu _ v )B (z) 



1 . i / u 2 - V 2 2K + ™)e~ 2 ™\ 



-(t; 3 - 3to 2 + 2(xv - uu x )) -3(ii(n 2 + -u 2 ) + 2(vu x + xu))e~ 2niu 
3(u(u 2 + v 2 ) + 2(w« a; + xn))e 2 ™ i> 3 - 3ra 2 + 2(xv - uu x ) 



+o(r 3 ) 



M(z) (8.9) 



as s — > oo and < r± < \z\ < r 2 < 1 (so \(\ — > oo). Since the function ((z) is of order 
O (s 1 ) on the latter annulus and Bq(z) is bounded, equation (|8.9|) yields the matching 
relation between the model functions U(z) and M(z), 

U(z) = (I + o{l))M(z), s^oo, < n < \z\ < r 2 < 1 

which is crucial for the successful implementation of the nonlinear steepest descent 
method as we shall see later on. This is the reason for chosing the left multiplier Bq(z) 
in (|8.7|) in the form 



9. Construction of a parametrix at the edge point z = +1 
Fix a small neighborhood U of the point z = +1 and observe that 
0(z) = 0(1) + i (4 + J) (z - 1) + 0{(z - l) 2 ), 
&s z £U and from (|6.6p 

S(z) = 0(\n(z- 1)), 

Both observations suggest to use the confluent hypergeometric function U (a, b; £) for 
our construction. This idea can be justified rigorously as follows. Recall that the listed 
confluent hypergeometric function is defined as unique solution to Kummer's equation 

zw" + (b — z)w — aw = 

satisfying the asymptotic condition as ( — > 00 and — ^ < arg £ < ^ (see [3] ) 

u(a , k , c, „ r » (1 _ !*I±J^ + -(. + D(l + ^ + + 0(( - 3) ) 

Secondly, using the notation U(a,C) = U(a, 1;C), the following monodromy relation 
holds on the entire universal covering of the punctured plane 

27T? 

17(1 - a, e-C) = e 2ma U(l - a, e"- C) - ^2(1 - a) ^ C)e ~ ? (9 - 1} 



and finally we have an expansion at the origin (compare to (|3.3p ) 

U(a,Q = c + c 1 ln( + c 2 ( + c 3 (ln( + O(( 2 ln(), C^O (9.2) 
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with coefficients Cj given as 
1 



co 



r(a) 



(</>(a) + 2 7E ), ci 



r(a 



(^(a + l) + 2 7s -2), c 3 = 



r(a) 



r(a) 



where 7^ is Euler's constant and = r(z)^ ■ Remembering the latter properties we 
now introduce the following matrix-valued function on the punctured plane £ S C\{0} 
(cf. [28]) 



Pch(C) 



U(u, e^C)e 27riiy e _i 2 



—tt < argC < 7T. 



xe 2V2 



(9.3) 



Let us collect the following asymptotic expansions. First in the sector — ? < arg £ < 



/ + 



( -4(i + ")' 



T(l-u) „ niu \ 

r(u) e 



c 2 V-W^( 1 + ^ e 



2„—iriv 



p(l-z,)V 

■4(i -v) 



+ o(r 3 ) 



For another sector, say \ < arg £ < we use ()9.ip in the first column of (|9.3p and 
obtain instead 

Pch(C) = 



+- 



-4(i + v) 



" r(i/) e 



c 2 V-¥F#d + -) 2e 



2„— nils 



4(1-^ 



+ o(C" 



xC 











1 



C -> 00, 



* 'J \ r(i-,)rw 7 
as well as for — =¥■ < arg C < — f with a similar argument in the second column of (|9.3 



Pch(C) 



1 



4 
,2 



r(i+^) -w 
> r(-i>) e 



r(i-i/) „ 7r ^> 



-1/ 



4(1 



#(1 + ")V 



r( 



4(1 -^) 2 



+ o(c 
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Also of interest is the following exact monodromy relation 



PcH(C) = PcH(e- 2ni C) 



g -2-KIV _|_ 



2tt 



27re 3 ' r " / 



(9.4) 



Let us now assemble the model function 



p§§ (0 = { 



e-^ fe 
e* 4 " 



^ch(C)I 

which solves the RHP depicted in Figure [6j 



(1-7)" 




arg C G (f ,tt), 
arg C E (-7r,-|), 



argCe(-f,f: 



l o 



27re 17rly 

r(i-i/)r(i/) 



(9.5) 



r(i-"-)r(^) 
1 

Figure 6. The model RHP near z = +1 which can be solved explicitly 
using confluent hypergeometric functions 

In more detail, Pqh(C) possesses the following analytic properties 

• p ch(0 is analytic for ( G C\{arg ( = -tt, -f , f } 

• The following jumps are valid, orienting the jump contours as shown in Figure [6] 

(^f(O), = (^f(C))„(l-7)- <T3 , argC = -vr 



r(i-i/)r(v) 



arg C 



7T 



(^#(C)) + = (i$f(C))_(J J , argC 

and in virtue of the classical identity 

r(i - u)r(u) = 

the entry of the latter triangular matrices equal 

27re~ i7TI/ 

i 7 (l -7)" 1 . 



TT 
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• As C — > oo, the model function Pch(0 snows t ne following asymptotic behavior 
in the whole neighborhood of infinity 



rgH(C) 



• / ..2 r(l-i/) 

i I v — WT 



r(i+v) -mv 2 



The model function P§§(C) will now be used in the construction of the parametrix near 
z = +1. Define 

C(z) = -2is 3 (tf(z) - 0(1)), \z-l\<r, (9.6) 

and notice 

C(z) = 2s 3 (4 + J) (z - 1) (1 + 0{z - 1)) , z -> 1 

yielding local conformality. This change of variables allows us to define the right 
parametrix V{z) near z = +1 as follows 

V(z) = ate-^Brizy^-^e-^W^Pgg (((^edCW+^W)^'!^^ (9 7) 

with ((z) in (|9.6p and the multiplier 

B r (z) = (c(z)^\y\ Br(l) = (l6s 3 + 4* S r : \ 

ijcr 3p -s 3 ^(z)a 3 



Also here, following from analyticity of B r (z) and conjugation with <j\e l i a3 e 

(1-7) CT3 



V p s 3 $(z)a 3 ( , 1 , 0L-A(zH 

X c ^ 7 (i_ 7 )-i i;c 

FIGURE 7. Transformation of parametrix jumps to original jumps 

in (|9.7p . parametrix jumps and jump curves (see Figure [7]) match original jumps and 
locally original jump curves. Furthermore, and we shall elaborate this in full detail very 
soon, the singular endpoint behavior of parametrix U{z) matches (|6.6p . i.e. 

V{z) = 0(ln(z-l)), z->+l. (9.8) 

Hence the ratio of S(z) with U(z) is locally analytic, i.e. 

S(z) = N r (z)V(z), \z~l\<\ 
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and again the role of the left multiplier B r (z) follows from the following asymptotical 
matchup 



V(z) 



if (i-«0<*3 p -S 3 tf(l)<X3 



v 



1 

xC 



-f(i + ^) 2 

-va 3 ^$(1)03 e -q(±-v)o 3 ^c 



r(l+z/) — ?rii^ 

, r(-*o e 

2 „iriv^ 



T(u) e 



+ o(r 3 ) 



+ 



4(i -vf 



r(i+^) ^2s 3 
r(-i>) 



+o(C 



M(z) 



"4(1 + -) 2 



(9.9) 



as s — > oo valid on the annulus < r\ < \z — 1| < r 2 < 1 (hence |£| — > oo) using the 
abbreviation 

If we are dealing with the case 7 < 1, then 

^(z) 1 - = O (s- 3±6Re ") =o(l), s^oo. 



This would mean that equation (|9.9p yields the matching relation between the model 
functions V(z) and M(z), 



1 



V(z) = (I + o(l))M(z), s-^00, < n < \z - 1| < r 2 < - 



(9.10) 



which is again crucial for the successful implementation of the nonlinear steepest descent 
method. However, if 7 > 1, then 



and hence 



/3 r 2 (z)i = /3 r 2 (z)^±j = 0(l), 



CO; P r (z) = C{z) 



Z + l 

z-1 



I'D 



With this notation, we have to replace (|9.10p in case 7 > 1 by 



V(z) = E r (z)(l + o(l))M(z), 



OO, < T\ < \z — 1| < T2 < 



where 



E r (z) = I _.rxw) e _ 2 ^ (1) ^ 



(9.11) 
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The appearance of the nontrivial matrix term E r (z) instead of the unit matrix in estimate 
(|9.1ip yields a very serious change in the further asymptotic analysis comparing with 
the matching case (|9.10p . We will proceed with this analysis in sections [13] and Q3J 



10. Construction of a parametrix at the edge point z = — 1 

-1. Opposed to 



For now, we introduce the model RHP near the other endpoint z 
consider 



Pch(() 



C/(-^e-^C)e"^ 



wiujj r(i+f)' 



xe 2 ^2 



<T2PcH(e- i7T ()cT2, < arg C < 2vr. 



and 



f ~ ( 1 o 

Pch(C)[ i 

Pch((){ 1 T(i-")rM 
V 1 



e l 7 





e-^" 



Pch(C) 





o ' 



arg C G (0, f ), 
argCe(^,27r), 



argCe (f ,f) 



The model function i^Jr (0 s °l ves the RHP of Figure [8] 



1 o 

27re-' ,r,y -, 



r(i-i/)r(i/) 
1 




(1-7)" CT3 



(10.1) 



Figure 8. The model RHP near z = —1 which can be solved explicitly 
using confluent hypergeometric functions 



• p ch (0 is analytic for ( G C\{arg C = 0, ^, 4f } 

• Along the contour in Figure [8l the following jumps are valid (recall (j9.4[) and the 
symmetry relation Pch(C) = a 2PcH (e~*' r C) cr 2) 

(^f(C)) + = (^f(C))_e- 2W3 , argC = 

(^#(0) + = (^#(0)_f_ J— f), arg C =^ 

\ r(i-i/jr(i/) 7 J 

(^(C)) + = (^f(C))_f? -^7TOV arg C = ^ 
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From symmetry P§§{Q = cr 2-f 5 Ci?( e_i7r C) fJ 2 and the asymptotic information de- 
rived earlier for Pch(C) m the different sectors, we deduce the following behavior, 
valid in a full neighborhood of infinity 



1 ( "4(1 -v? Tg?(l + ")V 



"c 2 -4(1 + -) 2 



+ o(C 



?(-!/) 

>V1 _ ,,^2^,711^ 

i» 

x(e- m C)" <T3 e- 4 l <T3 e^(5-)-3 ) 

Now, similarly to what we did in the construction of V(z), define 

C(z) = -2is 3 (${z) -tf(-l)), |z + l|<r (10.2) 

with 

C(z) = 2s 3 (4 + J) (z + 1) (1 + 0(z + 1)) , z -> -1, 

hence a locally conformal change of variables. With this change the parametrix W(z) 
near the left endpoint z = — 1 will be defined as 

W{z)=a 1 e^B l {z)e- i ^- v ^e- s3 ^ 1 ^ 

(10.3) 

with C(^) as in ([1021) and 

= \f~ W & z )j^l) > ^(-1) = (16s 3 + 4x S )" W3 . 

Since 

r(i-,)r» =Z7(1 - 7) 

the stated conjug ation with a^-^e-" 3 ^^ will again match parametrix jumps with 
original jumps locally on the original jump contour (see Figure[9j). Moreover the singular 



e 



s 3 #(z)a 3 ^ 1 7(l- 7 )- 1 ^ e -s 3 i9(z)a 3 

(1-7)° 




c V -7(1-7)" 1/ 

Figure 9. Transformation of parametrix jumps to original jumps 
endpoint behavior 

W(z) = 0{\n(z + l)), z^-1 (10.4) 
matches (|3.3p . leading to a locally analytic ratio of S'(.z) with W(z) 

S{z) = Ni(z)W{z), \ z + l\<\. 
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Also here the left multiplier B[(z) provides us with an asymptotical matchup 



W(z) 



1 + 



1 



^(l-^'C-i)^ 



C 2 \- T ^^ + ^- 2sH{ - 1) PT 2 ^) 



+o{c 



M(z) 



(10.5) 



as s — > oo valid on the annulus < n < |z + 1| < r2 < 1 (thus |£| — > oo) and we 
introduced the abbreviation 



z - 1 



z + 1 

Similar to the previous situation this implies on the annulus for 7 < 1 

W(z) = (I + o{l))M(z), s^oo 

whereas in case 7 > 1 

W(z)=E l (z)(l + o(l))M(z), s^oo 

with 



(10.6) 



Ei{z) 



1 -i^e^^pf(z)^ 



1 / ' ' " s ' V JV 'z + 1 



z - 1 



At this point we can use the model functions M(z),U(z),V(z) and W(z) to employ a 
further transformation. 



11. Third transformation of the RHP 
We put in this transformation 



R{z) = S(z) < 



(U(z))-\ 
{W(z))-\ 
{ (M(z))-\ 



\z - 1| < n, 
M < r 2 , 
|z + 1| < n, 

1,2 — 11 > ri, l-z + 11 > ri, |z| > r2, 



(11.1) 



where < r±,r2 < \ is fixed. With Co, r ,z denoting the clockwise oriented circles shown 
in Figure [lOl the ratio- function R(z) solves the following RHP 

• R{z) is analytic for z € C\{Cb,r,z U Ui=i 7«} 

• For the jumps, along the infinite branches ji 

R+(z) = R^{z)M{z)e sSi, ^G s e- s3 ^ z)(73 (M(z)y 1 , 
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Figure 10. The jump graph for the ratio-function R(z) 



with Gs denoting the corresponding jump matrices from ()6.5|) . On the clockwise 
orientied circles Co and C r< i, the jumps are described by the equation 

R+(z) = R-{z)U{z){M{z)Y\ zeC , 
R+{z) = R^{z)V(z)(M(z)y\ z G C r , 
R+(z) = R.(z)W(z)(M(z)y\ zGQ. 

• R{z) is analytic at z = ±1. This observation will follow directly from (|9.8[) and 
(|1U,4|) , which will be proven in section [TU] 

• In a neighborhood of infinity, we have R(z) — > I. 

We emphasize that, by construction, R(z) has no jumps inside of the circles Co )? -,z and 
across the line segment in between. In order to apply the Defit-Zhou nonlinear steepest 
descent method to the ratio-RHP, all its jump matrices have to be close to the unit 
matrix, as s -> oo, compare [18] . Hence it is now important to recall the previosuly 
stated behavior of the jump matrices as s — > oo: As mentioned before, due to the 
triangularity of Sjj, Sl and the sign-diagram of Re ^{z), the jump matrices corresponding 
to the infinite parts Ui=i 7« °f the -R-jump contour are in fact exponentially close to the 
unit matrix 

||M e s3 ^a 5e - S M>3 (M) -i_j|| L <C1 { e-^^l from Co; 



as s — > oo with constants q > whose values are not important. Also by virtue of (|8,9p . 
U(z)^M(z) s j approaches the unit matrix as s — > oo, 

||C/(M) _1 - J|| L 2 nioo(Co) < c A s~ l (11.3) 

with a constant C4 > 0. The jumps on C r j however have to be treated more carefully. 
As we already mentioned in sections [9] and [TUl estimates (j9.9j) and (|10.5p yield in case 
7 < 1 

HV(M)- 1 - I\\ L 2 nLac{Cr) < c 5 s-\ Hiy(M)- 1 - I\\ L2nLOC{Cl) < c 6 s- 3 (11.4) 
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as s — > oo. The estimations (jll.2p . (jll.3p and (jll.4p . which are uniform on any compact 
subset of the set (|1.16p 

{(7,x) G R 2 : -oo < 7 < 1, -00 < x < 00}, 

enable us to solve the ratio-RHP iteratively in that particular situation. 

12. Solution of the RHP for R(z) via iteration, 7 < 1 

Let us denote with Gr the jump matrix in the ratio- RHP and with Hr the underlying 
contour. The stated RHP for the function R{z) 

• R(z) is analytic for z G C\T,r 

• Along the contour depicted in Figure [10] 

R+(z) = R_(z)G R (z), z g E r . 

• As z — > 00, we have R(z) = I + 0(z _1 ) . 
is equivalent to the singular integral equation 

R-{z)=I+^-. [ R-(w)(G R (w)-l)-^— (12.1) 

ZTTl J W — Z_ 

and by the previous estimates (|11.2p . (|11.3p and ()11.4p . we have 

\\Gr ~ ^llL 2 nL°°(E fl ) < (12.2) 

uniformly on any compact subset of the set (|1.16p . By standard arguments (see [15]). 
we know that for sufficiently large s, the relevant integral operator is contracting and 
equation (|12.1|) can be solved iteratively in L 2 (T,r). Moreover, its unique solution satis- 
fies 

\\R--I\\lHs r )<cs-\ s^oo. (12.3) 

The latter information is all we need to compute the asymptotic expansion for the 
Fredholm determinant det(I — 7-Kcsin) hi case 7 < 1. Before we derive the relevant 
asymptotics let us first discuss the situation 7 > 1. In this case 

Hy(M)" 1 - /|| L2nLOO(Cr) 0, \\W(M)- l \\ L 2 nLOO{Cl) -» (12.4) 

and we need to employ further transformations. 

13. Fourth transformation of the RHP - undressing 

The presence of the multipliers E r {z) and Ei(z) in (|9,lip and (|10.6p requires further 
transformations leading to a singular or solitonic type of Riemann-Hilbert problem. 
Following [21j 110] . we will show how to deal with the singular structure. 
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A key observation for our next move is that the jump matrices G r (z) = V(z) (M(z)) 1 
and G[(z) = W(z)(M(z)^ 1 admit the following algebraic factorizations 



GJz) 



E r (z)G r (z)- I _.r^l e - 2 ^(i)^ (0)| ±i 1 



Gi{z) 



Ei(z)Gi(z) 



- E fef i ^ {1) f ± i(2--1) 



+ o(r 2 ) 



(13.1) 



-^e 2 ^(- 1 )/3 / 2 (,)(l-2,) N 



c 







+ o(r 2 ) 



(13.2) 



as s — > oo and < ri < \z =F 1| < r-i < 1. We observe that ||Gr r j — I\\ — > as s — > oo; in 
fact, since \C(z)\ > cs 3 on C r U C\, we have that 

\\Gr,l ~ I\\LZnL°°(C r> l) ^ C 7 S_3 5 s oo. 

Hence the natural idea is to pass from the function R(z) to the function P(z) defined 
by the equations 

( R(z)E r (z), \z-l\ < n , 
P(z) = I R(z)Ei(z), |z + l|<ri, (13.3) 
{ R(z), \zt 1| > r x . 

with < r\ < | chosen as in ([11. ip . By definition, the function P(z) solves the following 



RHP: 



P(z) is analytic for z G C\(C , r ,z U {±1} U Uf=i 7t) 
P+(z) = P_(z)G P (z), where 

G r ,i(z), z e C T) u 

G P {z) = l U(z)(M(z)y\ zeC , 

k M(z)e s3,? W^G5e- s3,? W CT3 (M(z)) _1 , z e 7i , i = 1, . . . , 8. 

P(z) has first-order poles at z = ±1. More precisely let P(z) = (P^^(z),P^ 2 ^(z)) 
with denoting the columns of the corresponding 2x2 matrix valued 

function. We obtain from (flTTO , (flTHZ]) and (TCP]) 



res 2=+ i 



res z= _i 



P«(; 



P (2) (1 ) 



r(i - 1/) 
r(i/) 



-2^(1)^2 (1) 



p( 2 )(.) = pW ( -l ) ( 2 i^=^e 2 ^- 1 )A 2 (-l) 



(13.4) 
(13.5) 



• As z — > oo, we have P(z) — > I. 
First the latter four properties determine P(z) uniquely. 

Proposition 4. The stated singular Riemann-Hilbert problem for P(z) has a unique 
solution. 
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Proof. The residue relations ()13.4[) . f)13. 5|) imply 

1 



P {+) (z)^_2p_ J, \z-l\<r; 



1 / 

where P^\z) are analytic at z = ±1. Hence one establishes detP(z) = 1 via Liouville 
theorem using the normalization at infinity und unimodularity of the jump matrices. 
From this and representation (113. 6p . the ratio of any two solutions P\(z) and Pi(z) of 
the given P-RHP, i.e. 

P 1 ( 2 )(P 2 ( 2 ))- 1 , 

is an entire function approaching identity at infinity; hence Pi = P2, showing uniqueness. 

□ 

Secondly, all jump matrices in the P-RHP approach the identity matrix as s — > 00; 
however P(z) has singularities at z = ±1 whose structure is described by the residue 
relations (|13.4[) and (j!3.5|) . This type of Riemann-Hilbert problem is a known one in 
the theory of integrable systems. The way to deal with such RHPs is to use a certain 
"dressing procedure" which reduces the problem to the one without the pole singularities. 

14. Fifth and final transformation of the RHP - dressing 
We put 

P(z) = {zI + B)Q(z)(^ JQ, (14.1) 

where B € C 2x2 is constant and see immediately that Q(z) solves the following RHP: 

• Q(z) is analytic for z £ C\(Cb,r,J U IjLi 7») 

• Q+(z) = Q-.(z)Gq(z), where 



and 



G Q (z) = [*J VjU(z){M{z)Y l / ^J, zGCo 



as well as 
Gq(z) = 



• Q(z) — > I, as z — > 00 

The Q-jump matrix Gq{z) is uniformly close to the unit matrix: therefore the Q-RHP 
admits direct asymptotic analysis, which will be performed after we determined the 
unknown matrix B. Using the conditions (|13.4p and (|13.5p 

res, =+1 P«(z) = (I + B)QW(l) = (I + B)QW(l){-p), 

res 2= _!P( 2 )(z) = (-I + B)Q^\-l) = (-I + B)Q^(-l)(-p), 
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SO 

B - (q(D Q ,«<-« (?) ) (- 1 ?) (0(D Q ,«-!) (?))"' (14.2) 

where p was introduced in (|13.6j) . Let us see for which values of s the latter matrix 
inverse is well-defined. 



15. Solution of the RHP for Q(z) via iteration, 7 > 1 

Since 

\\Gq - I\\v*nL°°(p R ) < ess" 1 , s -> 00 (15.1) 

we can solve the Q-RHP via iteration. This problem is equivalent to the singular integral 
equation 

Q_( Z ) = I+,L f Q_^ w )(G Q (w)-l)-^- (15.2) 
ATXl J w — z~ 

which can be solved via iteration in L 2 (X#), its unique solution satisfies 

WQ- -I\\i?(ii R ) < cs _1 , s-^00. (15.3) 
Combining the integral representation 

Q{z)=I+^- [q_( w )(Gq(w)-I)—, z£X R 

ZTTl J W — Z 

with (|15.ip and (|15.3p . we conclude 

Q(±l) = I + 0(s- 1 ), s^oo 

Hence the matrix inverse in the right hand side of (|14.2|) exists for all sufficiently large 
s lying outside of the zero set of the function 

1 -p 2 

which consists of the points {s n } defined by the equation 

80 1, . , , , , , r(l - v) 7T 

-s„ + 2xs n + - ln(7 - 1) In (l6s J + Axs) - arg v = - + nvr, n = 1, 2, . . . 

O 7T 1 (Z/"J 2 

and which will eventually form the zeros of the Fredholm determinant as written in 
Theorem 11.191 From now on, when dealing with the situation 7 > 1, we shall assume 
that s stays away from the small neighborhood of the points s n . We have now gathered 
enough information to prove Theorem Q] and [2j 
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16. ASYMPTOTICS OF In det(I — 7-FTcsin) - PRELIMINARY STEPS 

In order to prove the stated theorems we will use Proposition 11121 and [3j which in 
particular requires us to connect X(zLs) and X'(±s) to the solution of either the .R-RHP 
or the Q-RHP, see and P~15j) . To this end recall (JBSJ), ([TTTT]) and Figure O which 
implies for \z — 1| < r\ 



R(z)V(z)L~ 1 (z)e s3 ^ a3 =X(zs) 
On the other hand for \z + 1| < r\ 

R{z)W{z)L~\z)e sH ^ =X(zs) 



1 + 



2vrz 



1 + ^- 

2-ni 



-i V U + i 



-i iy U + i 



(16.1) 



(16.2) 



This shows that the required values of X(±s) and X'(±s) can be determined via com- 
parison in (|16,1|) and (|16.2|) once we know the local expansions of V(z), respectively 
VF(2;) at z = ±1. Our starting point is (j9.2|) 

di(C, v)e 2mv -d 2 ((, 1 - v)e™ V -^\ 



Pch(() 



-d 3 ((,i + ,) e ™^ 



r(-i/) 



-^(Ci-^^^ 

d 4 (C,-^) 



+ 0(C 2 lnC) 

(16.3) 



with (recall <^72\i ) 

di(C,v) = co(u) +ci(i/)^lnC + ^J, d 2 (C,v) = c H + c\ (u) ( In C - ^ 

and 

daCC") = -^i(C,^+^C2(^)+c 3 (^)(lnC + i^ 

as well as 

di{C,v) = ^d 2 ((,v) -i(c 2 (u) + c 3 (j/)^lnC - 

At this point use the changes of variables 

( = Q( z ) = -2is 3 (tf(z) - \z - 1| < n, A = zs 

and deduce from ()16.3p for — ^ < arg(A — s) < ^ using ()9.5p 



Pi(ln(A-s)) + (A-s)P 2 (ln(A-s)) + 0((A - s) 2 ln(A - s)) 







where the matrix functions -Pi (A) = (P 1 J (A)) and P 2 (A) = [P 2 (A)) can be determined 
from ()16.3p . For the remaining sectors — ir < arg (A — s) < — ? and ^ < arg (A — s) < tt 
we can derive similar expansions, they differ from the latter only by multiplication with 
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a triangular matrix. Now we combine the last expansion with (|9.7p and (jll.ip . so as 
A — > s in the sector — § < arg(A — s) < |, the left hand side of (|16,ip reads as 



R(z)V(z)L~ l {z)e 



1 ( rV 8 ^! 2 )^ 



R(z)a 1 e- i ^ 3 B r (z)e i ^-^ s e- sH ^P§§(C(z)) 



xe 1 4^(71 



fi(l)aie-^ CT3 J B r (l)e i 5^-^) CT 3 e - S ^(i)<T3 Pl ( ln ( A _ s )) 



e- iS r" 



o e *f " e ^ 3<Tl + (A " s) 



e *i(3-«')^ e -Mi)^ip L (i n ( A _ a )) 

xe i f^-^ CT3 e- s3,? « CT3 P 2 (ln(A - *)) 
+0((A-s) 2 ln(A-s)). 



i2 , (l)<7ie -i 5° s B r (l) + ^(lJffie-^B^l) 



e 4 2 







e**" 73 ^ + R(l) ai e- l * a3 B r (l) 



In the other sectors — ir < arg (A — s) < — ^ and ^ < arg (A — s) < tt we can derive 
similar identities, they differ from the latter one only by multiplication with a triangular 
matrix, see (|9.5p . The right hand side in (|16.ip implies in the sector — f < arg(A — s) < ^ 



R{z)V{z)L- l {z)e 



1 ( A^ 3 ^)^ 



X(zs) 



X{s) + (A - s)X'{s) + (9((A - sf 



7 /'-l A, A-s 
/+ 27nU lJH^ 



an - 2^(^n( s ) + *1200) ~ «) ai2 + ^?(*ii00 + *12(«)) ln(A - s) 

,a 2 i - 27i (*2i (s) + X22OO) ln(A - s) a 22 + £,(X 21 (s) + X 22 (s)) ln(A - s) 

+ (A- S ) 

x fAi-^^W + l^^MA-s) /8i2 + ^(^iW+^( 2 (s))MA-s) 

V/521 - 2^(^21 («) + ^W) ln(A - s) /3 22 + 2^(X 21 ( S ) + X 22 ( S )) ln(A - a) 
+0((A-s) 2 ln(A-s)), A^s 
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with constants otij depending on Xij(s) and fyj involving both Xij(s) and X-j(s). Com- 
paring now left hand side and right hand side in (|16.ip we deduce 



X u (s) 



R^iBril))- 1 

R(l){B r (l)Y 



11 



12 



Cq(1 - I/) 



r(i - „) 



7T 



+ci(u) \n(l6s 6 + 4xs) - i- 



X l2 {s) 



+ 



X 11 (s)+X 12 (s) = —e i ^ 
7 



R(l)(B r (l))~ 
+ci(-v) ( ln(16s 3 + 4xs) + 
R(l)(B r (l))- 
hence in particular 

'RWiBrQ.))- 1 
Furthermore 



cq(1 + v 



1» 

r(-i/) 



3 i?(i) 



S 3 1?(l) ' 



12 r(i/) 



*21 («) 



^(^(l))" 1 



21 



e i^ e s 3 m ( Co (_ u } + Cl (_^ / ln(16s 3 + 4xs) _ i 



R0.)(B r 0))- 1 ] e^e-^Hcoil-u) 



22 



r(i-r/) 

1» 



+ci(v) ln(16s 3 + 4xs) - i 



7T 



and 

^22 (a) 



e ifi/ e fl 3 »>(i) 



Cq(1 + y 



r(i + ^) 
r(-iz) 



fl(l)(B r (l))- 
+ci(-z^)^ln(16s 3 + 4xs) + 

12(1) (J3 r (l)) _1 ] e*i "e-- 3 *^) (c (z,) + Ci (i/) f ln(16s 3 + 4xs) + i£ 



which implies 



*2i(*) + ^ 22 ( S ) = HE e *W [fl(l)(B r (l)) _1 



S 3 1?(l) 



2lT(-l/) 



S 3 tf(l) 



22 Y[y) 
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Comparing after that terms of 0((A — s) ln(A — s)) we obtain 
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7 

1/03 3 + 



R'{1)-R{1) 



Is 1 



2 1 + 



■i.s- 



R{l)(B r {v))- 1 



R'(l)-R(l) 
(B r (l)y l 

s 3 tf(l) 



I/C73 3 + 



Is 2 



9 1 -I 2- 



(^(l))- 1 



e s 3 tf(l) 

n sr(-i/) 



-s 3 #(l) 



11 V2 



r -!/ 



12 sI » 

R(l)(B r (l))- X 



+ (8s 2 + 2s) 



i \ e 
ii/ 

12 V 2 



-s 3 tf(l) 



and 



^i(-) + ^(«) = ^*f| 

7 



is*' 



21 



R'(l)-R(l) 



2 1 + 



21 V2 



Is 2 



+ IV 



-s 3 tf(l) 



r 



22 sT ( U J 
R^iBril))' 1 



+ (8s 2 + 2x) 



i \ e 
iv 

22 V 2 



-S 3 1?(l) 



Although we derived the previous identities from a comparison in the sector — £ < 
arg(A — s) < ^, the same identities also hold in the other two sectors. There one uses 
the correct triangular matrices in (|9.5|) on the left hand side. Thus all previously derived 
identities follow in fact from a comparison in a full neighborhood of A = +s. Moving on 
to a neighborhood of A = — s, a completely similar analysis provides us with 



2-7T7 

X u (-s) + X 12 (-s) = — « 

7 



s 3 tf(-l) 



12 T(-y) 



fi(-l)(A(-l)) _1 



ii rd/) 



and 



2vri 



X 21 (- S )+X 22 (-s) = 



s 3 #(-l) 



22 r(-i/) 



^(-1)(^(-1)) _1 



e s 3 i?(-l)' 

21 T(v) 
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Also, comparing terms of 0((A + s) ln(A + s)), 



10 



*'n(s) + X' l2 {-s) = 2 -^e^ ( [(R'(-l) - Lt^) ( Bl (-1)) 



i» 2 



e -s 3 #{-l) 



R'(-l)-R(-l) 



2 1+ x 



Is' 



(fll(-l))' 



e s 3 tf(-l) 

ii sT(y) 



+ (8s 2 + 2x)|[ J R(-l)(i? / (-l)) 

-[«C-DW-i))- 1 l ia (| + 



-1 



S 3 tf(-1) 



-s 3 tf(-l) 



r(-i/) 



12 



and 



2-7ri „• w , 



W 3 3 + X3 



*£iH0 + = — ^( [(#(-1) - fl(-l) YlTA J 



R 



'(-D-^-D^^IW-i))" 1 



s r(-i/) 

+(8s 2 + 2x)^[R(-l)(B l (-l)) 



s 3 tf(-l) 



21 a r(i/) 



/ 3i . \ e 

\-lV) ; 

22 V 2 J T(-v 



s 3 i?(-l) 



22 



We finish this section by evaluating the resolvent kernel at A = ±s. Recall (|4.3|) and 
deduce 



Fi{±s) 
F{(±s) 



^-{X u (± S ) + X 12 (±s)), F 2 (±s) = J^-.(X 21 {±s) + X 22 {±s)) 



<Xi u (±s) + X' 12 (±s)), F! 2 (± S ) = x j^-{X' 21 (± S ) + X' 22 {±s)). 



2m 



Since 



R(s,s) = F{(s)F 2 (s) -Fi(s)F 1 (s) 
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we use the previously derived identities and obtain 
2m' 



R(s,s) 



1 



R' n (l)R 21 (l) -i4 (1) (16s 3 + 4xs) 



-2u 



a 2s 3 tf(l) 

sT 2 (-u) 



+ 



— 2s 3 $(l) 

# 12 (l)i? 22 (l) - # 22 (l)# 12 (l)] (16S 3 + 4xs) 2 " 6 gr2 
# n (l)i?(l) 22 " ^ 22 (l)i?ll(l) + # 12 (l)#2l(l) " # 21 (1)# 2 (1) 



i?H (1)^22 (1) - ^21(1)^12(1))^^^ 



1 



flll(l)i?22(l)-i?2l(l)-Rl 2 (l) 



Is 

8s 2 + 2:r 



T(u)T{-u)J- 
To simplify this expression, we have 

Proposition 5. R{z) is unimodular for any ijGl, i.e. &etR{z) = 1. 

Proof. Since detP/^(C) = 1, we have detU(z) = 1. Similarly detPg§(() = 1 = 
det Pqjj ((), leading to det V{z) = 1 = det W(z). Thus the ratio RHP has a unimodular 
jump matrix Gr(z), which means that the function det R(z) is entire. By normalization 
at infinity we end up with 

det R(z) = R n (z)R 22 (z) - R l2 {z)R 21 {z) = 1, z G C. 

□ 



In light of the last proposition, we obtain 



R(s,s) 



,2m 8s 2 + 2x 



-le 



+ 



2m 



ve 



3 + 



is-' 



2m 



7 T(u)T(-u) 7 sT(v)T(-v)l + ^ 7 sF(v)T(-v) 
R' n (l)R 22 (l) - # 22 (l)#i(l) + # 12 (l)#u(l) - #21(1)^12(1) 
2m ;„ 



+ e" 

7 



2s 3 tf(l) 

R' n (l)R 21 (l) - # 21 (l)#i(l) (16s 3 + Axs)-' 



sT 2 (-v) 



+ 



# 12 (1)#> 2 (1) - # 22 (1)# 2 (1) (16s 3 + 4xs) 



2v 



e -2s 3 i?(l) 



sY 2 (v) 



(16.4) 



Similarly, using Proposition [5] once more, 

, 2m 8s 2 + 2x 



ve 



3 + 



2m 



E(-s -s) = -ie lwu -— — ' I- 

1 ' ' 7 mn-v) 7 sT(y)T{-v)l + ^ 1S T(y)T{-v) 

#n(-l)# 2 (-l) - # 22 (-l)#i(-l) + # 12 (-l)#i(-l) - # 21 (-1)# 2 (-1) 



7 



— 2s 3 i9(l) 

R' u (-l)R 21 (-l) - # 21 (-l)#i(-l)l (16s 3 + Axs) 2u ' 



sT 2 (v) 



+ 



#i 2 (-l)# 2 (-l) - # 22 (-l)# 2 (-l) (16s 3 + 4xs) 



e 2s 3 i?(l) 



sY 2 {-v))' 



(16.5) 
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17. ASYMPTOTICS OF In det(7 — 7-FTcsin) - PROOF OF THEOREM Q] 

The stated asymptotics (|1.15p is a direct consequence of Proposition [3l Since the 
relevant estimations (|12.2p and ()12.3j) were uniform on any compact subset of the set 
(|1.16p . the asymptotic series for 



d_ 



In det (I - Tocsin) 



can be integrated with respect to 7, which leads to (|1.15p including the constant term 
X- We trace back the transformations 

X(X) 1 y T(z) 1 y S{z) 1 y R(z) 

and deduce first asymptotic expansions for the coefficients X\ , X2 and X 3 in (|4.15p . First 



X l = lim (A X(A)e~ H 3 

A— >oo 



i{\\ 3 +x\) 



cr-.i 



Ij) = -2ua 3 s + ^ / R_(w)(G R (w) - I)dw 



-2va 3 s + 



is 
2^ 



is 



[R-(w) - l){G R (w) - I)dw + — I (G R (w)-l)dw (17.1) 



and we have for z £T, R from (fT2TTJ) . (fl2~3]) and (USD 



C 



ss/w-))- 1 



-ue 



v ue 

2niv 



Bq(w) 



idw 



Co 



2sw(w — Z-) 



1 

2s~z 



V u 
-u —v 



{Bo(z))- 1 



-ue 



v ue 



B (z) 



+ o(,- 2 ) 



+ o(,- 2 ). 



The latter expansion will be improved via iteration 

R-(z)-I = ±-J (R_( w )-l)(G R (w)-l) 



dw 



w — Z- 



C 



+-L f(G R (w)-l)-^-+0{8~ 3 ), s^oo. 
2in J w — z- v 



Co 

If [•, •] denotes the usual matrix commutator on C 2x2 , (|8.8p implies for any constant 
matrix A G C 2x2 

(B (w))~ 1 AB (w) = (B (0)y 1 AB (0) + 2vw[(B (0)y 1 AB (0),a 3 

+4u 2 w 2 \(B {0))~ 1 AB {0),a 3 \a 3 + 0(w 3 ), w -> 0. 
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This in turn leads to 



R-{z) - I 



2sz 



V u 
-u —v 



1 



+ 



4s 2 z 2 
1 



8s 2 z 2 



V u 
-u —v 



u 2 — v 2 



V u 
-u —v 



{Bo(z))' 1 
- (Bo(z))- 1 



-ue 



v ue 

2niu 



-2niv 



B (z) 



-ue 



v ue 

2niv 



-2iviu 



B (z) 



2(u x + uv 
2(u x + uv) u 2 — v 2 

u 2 -v 2 2(u x + uv)\ 
2{u x + uv) u 2 -v 2 )^ {Z) 

for any z 6 Back to (j!7.1|) . one starts with 

2iriv 



[Bo(z)y 



v (u 2 


-u x \ 




s 2 z \u x 


-u 2 ) 



+ 0(i 



oo 



h 



[R-(w)-l)(G R (w)-l)dw 



—u 

uv u 



uv 

2 



+ 



i{-2m)v 1 



2uu^ 



2u 3 - v(u x + uv) - tt(u 2 - v 2 ) 



-2u 3 + v(u x + uv) + §(tt 2 
— 2uu T 



+0(s 4 ), s — > oo 
and moves on to 



h 



Si, 

+ 
+ 



(G R (w)-l)dw 
%{-2-ki)v 2 



(— 2iri)i ( v 



2s 



-u —v 



+ 



-2*Kl)V 





u x + uv 



-{u x + uv) 




-(§(" + v ) + vu x + xu) 

|(u 2 + v 2 ) + vu x + xu 



(-2iri)i 
4s 3 + xs 



- V— v 2 cos a 



V 



2 ) 



with 



. O o 

cr = a(s,x,j) = -s +2xs + 



ln|l-7| 



7T 



In (16s 3 + 4xs) — arg 



r(i - v) 

1» 



(17.2) 



Together 



I1+/2 



(— 2wi)i ( v 



2s 



u 



+ - 



-2iri)iv 2 (2uu 



-u —v 



-2uu 



+ 



+ 



{—2m)v ( u 



s- 
-2iri)i 



u, 



4s 3 + 



xs 



cos a 



V 



cos a 

.2 



+ o(^- 4 ). 



where we used that it = u(x, 7) solves the second Painleve equation. The previous line 
leads us to the following expansion for X\ 

,2 .. \ „-,,2 



Xi = — 2z/a"3S + 
+ 



4s 2 + 



V 







2 






V- 


-Z/ 2 COS (7 





+ 



v lu 



s \u x 



u x \ iu f 2uu x 
2]+ ^\u xx 



-u 



-2UUn 



■v 2 co: 
^ 2 



+ 0(s 3 ), s ->• 00. 



(17.3) 
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Secondly 

,2 -2 



z ■ 2 r 

X 2 = 2u 2 s 2 I- l -^-(h+I 2 )a 3 + ^- R_{w)(G R (w) - I)wdw. (17.4) 



Si? 



We need to compute 

Is = J (R-{w)-l){G R (w)-l)wdw 



{-2m)iv ( u 2 v + uu x §(u 2 -t; 2 )\ 

2s 3 \u {u 2_ v 2 ) U \ + UU J +U \ S ) 



as well as 



r x , ( — 27rz) / u 2 - v 2 2(u x + uv 

(— 2-Ki)iv ( u(u 2 + u 2 ) + 2(una; + xu) 

4^3 \u(u 2 + v 2 ) + 2(vu x + xu) 



2iri ( \J — v 2 si 



TlF— n +°(* ) 

4s cl + xs l-y-f smff L) / v ' 



to obtain 



(— 2m) ( u 2 - v 2 2(u x + uv) 
-8^{2(u x + uv) u 2 -v 2 

(—2iri)ii> I u 2 v + uu x u 3 + vu x + xu 



2s 3 V ^ 3 + vu x + xu u 2 t> + uu x 



2-ki ( \J —v 2 sin (A , _/ _ 4 n 

4s ,r> + xs \V — ^ sin cr U / 



Adding together in (|17.4p 

,2„2 r „•„„/' « _M ^-2^ 2 f u2 U;r ^ + - f u2 ~' y2 2(li x +W) 
-u v ) \u x u 2 ) 8 \2(«a; + ««) u 2 — v 2 



2iu 3 ( 2uu x xu + 2u 3 \ iv ( u 2 v + uu x u 3 + hm^ + xu 



•S 



xu + 2u 3 2uu x / 2s V u d + fu x + xu u z w + uu a 



+ ■ 



2ivs ( v 2 y/—u 2 coscA s ( \J —v 2 sin (j 



4s 2 + x \\/— z^ 2 cos a z/ 2 y 4s 2 + x V yf—v 2 sin cr 

+0(s" 2 ), s^oo. (17.5) 
Finally the computation of X3 

X 3 = -^ s s (l + 2v 2 )a 3 + —{l 1+ I 2 )-—{l 3 +h)a 3 

■ 3 r 

/ R-(w)(G R (w) - I)w 2 dw. (17.6) 
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Since 



h= / (R_(w) - l)(G R (w) - I)w 2 dw = 0{s~ 4 ) 



E 



R 



and 



t-, i \ t s~i i -v T \ 2 7 (-27^^)^ / -zv 2 V— ^ coscA 

J 6 = R-(w)(G R (w) - l)w 2 dw = — , 5- 2 

{-2-ni)i ( -(v 3 - 3vu 2 + 2(xu - to x )) -3(u(ti 2 + v 2 ) + 2(to x + ra)) 



^ 48s 3 \ v 3(u( , u 2 + w 2 ) + 2(vu x + xu)) v 3 - 3vu 2 + 2(xv - uu x ) 

we add together and obtain in ()17.6p . using the previous results 

o 2 . 2 2 ( v u\ „ q /u 2 — 

X 3 = - yS 3 (l + 2,>3 + z,V^_ n _J+2, 3 S (^ _ n I 

z^s ^ n 2 — f 2 — 2(u x + to)^ ( n . ^(2uu x —u xx 

"x 



4 ^2(1^ + to) -(u 2 -t; 2 ) J+ 2?t/ y Uxx -2uu 
2iv 2 s 2 ( -v 2 \/-u 2 cos a\ is 2 ( -v 2 \f- 



^ 4s 2 + x \— V— v 2 cos cr z/ 2 / ^ 4s 2 + x ^— V 

2vs 2 ( —yj—v 2 sincA . 9/ ii 2 t> + to x — (u 3 + vu x + xu 



-v cos a 



' 2 , 
,2 



2 / (A f T" wu x — V u ~r v "a: 
,3 



4s 2 + x \\/-v 2 sin a / \« + vu x + — (irv + uu x ) 

i f-(v 3 - 3vu 2 + 2(xv - uu x )) -3(u(u 2 + v 2 ) + 2(vu x + xu))\ n , _^ 



48 \3(u(u 2 + v 2 ) + 2(««a; + xu)) v 3 - 3vu 2 + 2(xi; - uu x 
With the given information at hand, fj4.9|) and ()4.10p lead to 



2v 2iv 2 u xx 2i\f 



2 



b = 2X\ 2 = iu u x + — 1 cos a + (s" 3 ) (17.7) 

s s 4s^ + x 

as well as 

c = 2X 21 = -b + 0(s~ 3 ), s^oo (17.8) 

and 

d = ix + SiX{ 2 X 21 (17.9) 



2 8z/ 8zz/ , 9 2 4\ 8tu\/—u 2 , 

= «x + 2«n to x ^-((^x) + xu + 2u H — — = cosct + Gms 

s 4s 2 + x 

Furthermore 



e 



^(XPXf - X 2 12 ) = -2i« x + — {xu + 2u 3 ) + — sina + 0(s" 2 ) 

\ J S 4$ I 37 

and 

/ = _ e + ( s -2) (17 . 10) 

where we made use of the following identities, see (|17.3p and (|17.5p 
X 21 = -Xl 2 +0{ S ~ 3 ), Xl^-Xf+Ois- 3 ), X 21 = X?+0{s- 2 ), s^oo. (17.11) 
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We have now derived enough information to evaluate the first terms listed in Proposition 

El Put 

2 ^ ^ 



Vi{s,x,i) = -Si[{X z ) 1 + {X 1 ) 1 {X(-X 2 )-{X 2 ) 1 X 1 'j -2d[{X 1 
+4i&((X 2 ) 7 - {X^xX 1 - Uc( (X 2 ) 7 - (X^Xi L 



and notice that 



From ()f 7,f Op and (|17.f ip we obtain therefore 



-eUx^Y 1 = -f((X 1 ) y ) 12 + 0(s- 2 ) 



and via (fTTTfl) and (117.111) 

4i&((X 2 ) 7 ) 21 = -4ic((X 2 ) 7 ) 12 + O . 



Since also 

21 



-tibUX,)^)- -4ib((X?) y X? + {xf)xf 



4*c((xf) Xf + (X^) 7 Xj 3 )+o(^) 



as s -> oo uniformly on any compact subset of the set (|1 . 16|) . we can simplify the 
expression for Pl(s, x,j) asymptotically 

Vi{a,x,i) = -8i((X 3 ) 7 +(X 1 ) 7 (X 1 2 -X 2 )-(X 2 ) 7 X 1 ) 11 -2d((X 1 ) 7 ) 11 

+8z6((X 2 ) 7 - (X,)^) 21 - 2e((X 1 ) 7 ) 21 + o(^fj . (17.12) 

Next from (fT73I) 



X? 



22 2 2 n — v 8iv Aiv^s iv 2 

4u s — 2iuvs — 4u u -\ uu x + — ^ 1 [vu + uu x 

4 s 4s z + x s 
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with I denoting the 2x2 identity matrix. Thus 



17 



11 



+o 



In s 
~^3" 



Av 2 s 2 - 2ivvs - 4z/ 2 u 2 + 



2 (utl^^) 

-\ + — + -4 (2v 2 uu x ) ^ (v 2 ) 

2 7 s s 2K >i 4s 2 + x v ; ^ 

4zi/ 3 s 



2v^s + -Uy + 



u 2 — v 2 8iv 3 



-uu x + 



4s 2 + x 



IV 



+—{vu 2 + uu x ) + 0(s 2 ) 
s 



iv 2 v 7 s 3 + Aivv^vs 2 + 2iv 2 s 2 v 1 + 8v 2 v 7 u 2 s — v^s 



u 2 — v 2 



+ vsvv 7 



+4v 2 s(vu 2 ) + 16^^ 



uu q 



8ZZ^ UryS / 9 \ 9 9 

- 2 2ivv 1 [yu + uti x J — 2zf u v 7 



i u 2 — v 2 



+ 2 V "— 



2ivv[vu 2 ) + Aiv 2 (2v 2 uu x ) 



Moving on, we use 

{[x l ) i x 2 + {x 2 ) i x 1 

and obtain 

((x 1 ) 7 x 2 + (x 2 ) 7 x a 



11 



Mv 2 s 2 
i 4s 2 + x 



12 \ v 21 



2 ) +0 

1 1 



In s 

s 



-12\ v 21 



11 



4i/ 3 s 3 + 3iv 2 s 2 v + 6v 3 u 2 s - 3iv 2 u 2 v + 12iv 4 uu x 



6iv 4 s 2 



— IV uu,, 



vs- 



4s 2 + x 

2 i u x + uv 
-IV u x u 7 + -ii 7 



u 2 — 3v 2 iv u 2 — v 2 
-. h — 



vs s , N 

—uu 7 + -u[vu) 



ivu[vu x ) +ivu x \yu) + iu[y u x ) — — u 



'7 2 



+0( S " 1 ). 
Furthermore 









r 




— ivs 1 









AvVjS I — iVjS 



1 /t!tty — UU^y UV^y — VUry 



-uu v. 



2vv^ 



-U r U 



U x 
2 



4 \ uv -y ~ vv "y vv 'y ~ uu 



2v 



(VU 2 ) (^)7) +0(s -l) 

[yu x ) 7 (yu') 7 J 



and thus 



( \ 21 

( [Xij — \X-ij X\\ = is(vu)j + ivjsu — ivsu 7 — 2{v u x ) 7 + 



(u x + uv)- 



-2vv 7 u x + -(mv 7 — vu 7 ) + 2v{vu x )~ 1 + 0(s x ) . 



which implies 

b([x 2 )^-{x 1 ) i x 1 



21 



-suyvu) 



su 2 v 1 + svuu 7 — 2iu(y 2 u x ) 1 + iu 



(u x + uv) 7 



UVj — vu 



—2ivv~,uu T + iu- 
+2iv 2 u x u.y + O^s" 1 ). 



7 



+ 2ivu(vu x ) 7 — 2ivu x {vu)~ i — 2ivv 7 uu x 
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Finally 



and 



Let us write 



11 



-2iu.ysx — AiujU s — - — u u 7 + \Qvv^uu x + U[s ) 



21 



+ 0(s-i). 



V x (s, x, 7 ) = s 3 r[ 3) (x, 7 ) + s 2 r[ 2) (x, 7) + sV[ 1] (x, 7) + v[ 0) (x, 7) + O 
where P-^ (21,7) are independent of s. Since 



In s 

s 



17=0 







we get from (|17.12p and the previous computations 



7 7 r 

J V[ 3) (x,t)dt = J -8^-y(l + 2^ 2 )J -8i(-8v 2 v t ) +8i(-iv 3 ) t 



Next 



" r 

/ 




L 




16. 




IV. 


3 






T r 


J 
c 


- 



8i(zi/ 2 w) t - 8i(4ivv t v + 2i^ 2 u t ) + 8i(3iz^ 2 t>) t 



r/f 



0. 



and 



8i[ 2t 3 u 2 - -(u 2 - v 2 )j - 8i\8v 2 v t u 2 - v t 



+vvv t + 4v 2 (vu 2 ) t ^j + 8i^6v 3 u 2 - v 



u 2 - 3v< 



V u 
f 8i[ - -uu t + -[yu)t 



+8i \ iu(i(vu)t + ivtu — ivut)\ — 2 f — 2ivtx — AivfU' 
= Aivx. 

We conclude the computations for "Pi (s, 2,7) by evaluating 
7 7 - 



dt 



J v[°\x,t)dt 



4 2iu A S 2 is 2 V 2 0/2 \ 

IV (u V + uu x ) 
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(v 3 — 3vu 2 + 2(xv — uu x )) I — 8i( 16iv 3 v t uu 



4s 2 + x 4s 2 + x 

8iv 3 v t s 2 



4s 2 + , 



2ivv t (vu 2 + uu x 
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99 . U V , 0\ n • 2 / 2 \ 4zZ^ S / 9 \ 

-2?z/ u v t + iv t 2ivv[yu ) t + 8iv (y uu x ) t - — ^ \v \ 

8 <is z + x 

2 2 in- 4 -2 6iv 4 S 2 . . U 2 — V 2 



+8i\ — 2>iv u v + 12iz/ uk t — ii/ to t h ii> 

4s 2 + x 16 

+8i{ — iv u x Ut + iut ivu[yu x )t + ivu x [yu)t + ^^^(^ / u^Jt 



. (U x + Uv) t \ ( XV t 2 ,!« \ e Y o- / 2 N 

— zu 1—21 u v t + Vovvtuu x I + 8z I — liu[y u x )t 

(u x + Ut> )t Ut>t — TO* 

+iu 2ivvtuu x + iu h 2ivu(yu x )t — 2ivu x (yu)t 

—2ivvt uu x + 2iv 2 u x u-t^j — 2^ — ti^ii^ (it 

2 2 Z" 7 
= —2v H — — uMj;) + 2 / u x utdt. 
3 Jo 

The following observation is useful 

Proposition 6. Let u = u(x, 7), 7 < 1 denote the Ablowitz-Segur solution of the second 
Painleve equation u xx = xu + 2ti 3 as described in Theorem [ZJ Put 



7 

7) = ^(^(x,7) - u(x,j)u x (x,j)) + 2 y u x (x,t)u t (x,t)dt, 



F(x, 7) = - / (2/ - a;)w 2 (y, 7 )dy. 



F{x 

then 



Proof. Using the differential equation for u and integration by parts, we obtain 

—F(x,j) = v{x,i) = (u x {x,i)) - xu 2 (x,j) - n 4 (x,7) 
ox ' 

and hence after integration 



F(x, 7) = - J (y - x)u 2 (y, 7 )dy + C( 7 ) (17.13) 

X 

with a constant C, only depending on 7. Since u decays exponentially fast as x — > +00, 
the same limit on both sides of ()17. 13|) gives us the stated identity. □ 

At this point we summarize the previous computations. As s — > 00 uniformly on any 
compact subset of ()1.16|) 

7 00 

J Vi(s,x,t)dt = wj^ys 3 + isx^j - J (y-x)u 2 (y,j)dy + 2(iv) 2 + o(^^J (17.14) 
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To move further ahead let us define 

V 2 (s,x,j) = (A 11 -A 22 )H 11 (s) + A 12 H 21 (s) + A 21 H 12 (s), H(s) = — (s){X(s)y 1 
with (compare (|4.8p ) 

Since 

An = + ^W) (^2l(s) + X 22 ( S )) 

we can now use the identities derived in section [16] for Xij(±s). With 

m±1) = 1 + h I r - [w){Gr(w) ~ ^vri (17 - 15> 

. /+ _L/ (GBH _ 7) J|_ + 0(s - V/±i (_ : ») +0( .->) 

Co 

and the classical identity 



r(z)r(l - z) = — ^— , z€C\5 
sm7rz 



one concludes 

An = i/ + 0(a _1 ), A 22 = -v + 0{s- 1 ) 
uniformly on any compact subset of the set (|1.16p . Also 

A 12 = v(16s 3 +4xs)- 2 »e 2s:i ^^^ + 0{ S - 1 ) (17.16) 

and 

A 21 = -u(16s 3 + ^a) 2 ^- 2 ' 3 ^ 1 )^^ + Ois- 1 ). (17.17) 
Next we use (|17.15p to simplify the identities for Xij(s) obtained in section [TBI 
Xn(s) = (16s 3 + 4xs)-V5^( 1 )( Co (-i/) + Cl (-zy)^ln(16s 3 + 4xs)-i^ 



W— 



X 12 (s) = -(16s 3 + 4xs)- v e i ^e sH{l HcQ{l + v) T( > l / + V } +c 1 (-u)[ ln(16s 3 + 4xs) 

T(-v) 



+*-) 



X 21 (s) = -(16a 3 + 4xs)V^ e - s3,?(1) (c (l - " f^,^ + ci(^) ln(16g 3 + 4xs) 



. 7T \ \ / In s 

- ! 2 ) +0 (— 
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\-c l {v)( ln(16s 3 + Axs) +i| 



X 22 (s) = (16s 3 + 4xs) Vf V 
Combined with (|17.16|) and (|17.17p . we deduce the following asymptotics for V 2 (s,x,'j) 



T 2 (s, x, 7) = 2ue inu (co(y) + ci (1/) f ln(16s 3 + Axs) + 

c (-i/) + ci(-z/)^ln(16s 3 + 4xs) + (co(-f) + ci(-i/) f ln(16s 3 + 4xs) 



- 2zV™ (co(l - ^)^rp + c i (") ( ln ( 16s3 + ~ «|) 
ln(16s 3 + 4xs)^ (c (l + »/) ^ + ci(-i/) (Sn(16s 3 + Axs) + i^Q 



7T 



■ r(iz) 



+ ( c °( 1 + v) T T{-u) + c i("^) ( ln ( 16s3 + 4xs ) + *|) 

cq(i/) + ci(i/) f ln(16s 3 + 4xs) + ) (c (l - ^ ^7^ + c i0) ( ln(16s 3 + Axs) 



IT 

-«2 I 



1» 



x c (l - f) 



r(-i/ 
r(i - 



co(^) + ci(v) ( ln(16s 3 + Axs) + i^- 



i» 



+ ci(i/)^ln(16s 3 + 4xs) -*'^J 



-r(-v) 



i» 

( ln(16s 3 + 4xs) + 



7T 



co(-i/) + ci(-i/) ln(16s 3 + Axs) - i- I c (l + 1/ 



r(i + 

r(-iz) 



+ ci(-i/ 



7T 



co(-^) + ci(-i^) ln(16s 3 + 4xs) - i 



7T 



x (^C 1 + ^ri-J) + ^("^ ( ln (16s 3 + Axs) + i|J) + o 
What is left in the identitiy stated in Proposition [3] is the term 



(Ins) 



s —> 00. 



with 



and 



Vz{s,x,i) = {B ll -B 22 )H ll {-s) + B l2 H 2l {-s) + B 2l H 12 {-s) 



B u = u + Ois" 1 ) , B 12 = u(l6s 3 + Axs) 2 »e- 2sH ^^^- + C^s" 1 ) 
B 21 = -v(Ws 3 + Axs)- 2 ^ 2 ^^^ + Ois- 1 ), s -> 00 



r(-rz) 

which, also here, holds uniformly on any compact subset of (|1.16|) . Again (j!7. 15[) allows 
us to simplify the idenitites for Xij(-s) obtained in section [TBI and we are lead to the 
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V 3 (s, x, 7) = 2ve i * v (co(-u) + ci(-i/) f ln(16s 3 + 4xs) - 
+ ln(16s 3 + 4xs)^j (cq(v) + ci(i/) ^ ln(16s 3 + 4xs) + *^J) + ( c o W + ci( 



x ^ln(16s d + Axs) +i~ 
x ( ln(16s 3 + Axs) + 



2^(co(l + i/)^-t^ + Cl (-y) 
; i| + ln(16s 3 + 4xs)^) (c (l - I/ ) r( r(~) t/) 



+ Cl (^ln(16s 3 +4xs)-i^) + (co(l-i/) r( ^^ 



+ci(i/) ln(16s 3 + 4xs) - i- 



1» 



Cq{-v) + Cl(-i/) 



x ^ln(16s 3 + Axs) - ) (c (l + + <*(-*/) ^ln(16s 3 + 4xs) + 

W^^=^(co(-^) + Cl (-z,)(ln(16 S 3 + 4x S ) -z|^( Co (l + 
+ci(-z>)^ln(16s 3 +4xs) +«|J) + ^e^^TT ( c o0) + c iH f ln(16s 3 + Axs) 



. 7T 



+^0 Mi-" 



r(i-«/) 
r(i/) 



r(-rz) 



+ ci(i/)( ln(16s 3 + 4xs) - i 



7T 



^™r7=S) ( Co ^) + Cl ^) ( ln ( 16sS + 4xs ) + *f ) 

(cod - + dM (ln(16 S 3 + 4x S ) - ,|)) + 



, s — > 00. 



The two expansions for T > 2(s,x,-f) and ^(s, 33,7) combined together allow us now to 
evaluate 

7 



J (V 2 (s,x,t)-V 3 (s,x,t))dt. 



(17.18) 



For this evaluation it is important to recall the definitions of cq(v) and ci{v) 



1 



1 



as well as the functional equation of the Digamma function (see e.g. [3]) 

ip(z) = %/,(z + l)-- = </>(! - z) - vrcot vrz, z G C\{0, -1, -2, . . .}. 
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It implies 

cq(v)ci(-v) + a{v)co{-v) - co(l + v ) V y{-v) Ci ^ ~ Cl (~ l/ ) c °( 1 ~ v ) ~T(z/) = ° 

and shows therefore that all terms of 0((lns) 2 ) in (|17. 18f) vanish. The remaining terms 
of O(lns) and O(l) can be computed in a similar way, we obtain 

(V 2 (s,x,t) -T> 3 (s,x,tj)dt = 6(iz^) 2 lns + 8(i^) 2 ln2 (17.19) 

+2 J u(t)(lnT(u(t))-lnT(-u(t))) t dt + o(^^^ 
o 

as s — > oo uniformly on any compact subset of the set (|1 . 16)) . The latter statement 
combined with (|17.14|) implies (|1.15p with an error term of 

In order to improve this error term, we use Proposition [T] and in particular the identities 
for i?(±s,±s) derived in section flUl With (|17.15j) after simplification 

R(s, s) = -iv(8s 2 + 2x) - + 0(s~ 2 ) 



and also 



which implies via 



R(-s, -s) = -iu{Ss 2 + 2x) - + 0{s~ 2 ) 



\2 



^-lndetf / - -fK csin ) = iu(16s 3 + 4x) + 6 ^ + 0(s~ 2 ), s -)• oo 
as s 

uniformly on any compact subset of the set (|1.16|) . Integrating the latter equation with 
respect to s and comparing with (|17.19|) , we complete the proof of Theorem [TJ 



18. ASYMPTOTICS OF lndet(J — 7-ffcin) " PROOF OF THEOREM [2] 

The statement on the asymptotic distribution of the zeros of det(/ — 7-K" C sin)>7 > 1 
will follow from an asymptotic expansion for the latter quantity which we derive from 
Proposition [1] and [2j To this end we trace back the transformations 

X(X) t-> TO) t-> S(z) i — y R(z) i — y P{z) t-> Q(z) 

and use the identities (|16.4p . (|16.5p . which were derived independently of the choice of 
7 £ M. First the residue conditions (|13.4p and (|13.5p show explicitly that R{z) is analytic 
at z = 1 

*(l) = (g(l) + (I + fl)tf(l))Q °^+{I + B)Q(l) LJ_^L i) (18-1) 
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where vq and p were introduced earlier as 

"° = h I 1 - * = Si lnh " ^ ' = 

Furthermore 

rf(i) = WD + ^+WdW^li J) + (o'd) + + J) 

+<™< i >Ui, :C ) A) (i8 - 2> 



where 



Similarly 



x) = \(f (i + + (-0 - i) (3 + 2 ) (l + ^ 



^_fr_ R W_i^4 -^Pt^ 



-(I-B)Q(-1)I ^ ^T+^l (18 . 3) 



v 1 

and 

3+ 



R'(-l) = (q(-1) - (/- fl)Q'(-l)) ( ^ ^J+Sj (18.4) 

+ (q'(-D-(/-5)^)(2 j)-(i-B)0(-i)(7 ^f'^ 



To evaluate Q(±l), we iterate. First for any z 6 from (|15.1|) . (|15.3|) and residue 
theorem 



= i + ^Jq-H{g q (w)-i) 



dw 



w — Z- 







1 + 




2sz I V 









where we used that 
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Thus 



« ±1 ) = /+ i/(«-W- J )( G «W-/)^ + i/(G (w)-/)^ 



I± 



V ~U\ UQ -vr -u x 

¥ ~2 i U r U 2 



u 2 — v 2 



1 



8s 2 \-2(u x + to) 



2s \u —v j 
+0(s -3 ), s —?■ oo, 
uniformly on any compact subset of the set (|1.18j) , and similarly 



-2{u x + uv] 
u 2 - v 2 



Q'(±i) 



Since 



2s V u 
+0(s~ 3 ) 



v —u\ I —u 



-u x 

,,2 



1 



u 2 — v 2 



T 4s 2 V -2{u x + to) 



le 



with cr = a(s,x, j) as in (j!7.2|) . we obtain for the matrix 



p = p 



-2(u x + to) 
u 2 — v 2 



)( 1 )• 



which appears in (|14.2p . that 



v u 

det = — 2ip( cos a sin cM — 

s s 

H — [u x + u sm a J + 



u 2 — u 2 



2s2 cos ( 7 + 0(s- 3 )), (18.5) 

Since we agreeded that s stays away from the small neigbhorhood of the points {s n } 
defined by 

cos <r(s n , x, 7) = (18.6) 
we see that for all sufficiently large s lying outside of the zero set of the latter transcen- 
dental equation, the stated determinant is non zero. Back to (|14.2|) . this implies 



B 



2p 
dete 



sin a —1 \ coscr 
1 — sin a J s 



H — 2 { 2iv§ cos a 



-u —t 

u T u 



1 

+ 2 



u 

u —v 

—v 2 sine 
u 2 

s — > 00 



v sm a 



— 1 sin a 
sin <t 1 



and with 
2p 



detG 



cos cr 



1 H — ( f tan a 



COS (7 



H — n\ v 2 tan 2 cr 



2to sin a u 
+ 



COS 2 (7 



cos 2 a 



u 2 — v 2 



2iu u 2 tan a - 2lu ° Ux ) + O (s -3 ) 

COS (7 
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we obtain in turn 



B 



1 



cos a 



sincr —1 
1 — sin a 



1 f v sin a — u /sin cr 
1 



cos^ a 



sin cr 



+ 



v —u 
u —v 



1 



+- 



1 



2cos<7 

2 



2v 2 sin a — u 2 sin cr — u x 

2 



u — u x sin cr 

2iz/ (« 2 sincr + li^.) /sincr —1 

— sin <7 



cos 2 cr 



1 



+ 



—v + u x sm<7 
— 2v 2 sin a + u 2 sin cr + u x 

(v sincr — u) 2 / sincr 



cos 3 a 



1 



+ 



-uv cos cr — 2ll>QU 2 U 2 COS O" — 21uqU x 
-u 2 cos cr + 2iuQU x uv cos cr + 2iv§u 2 



+ 0{s 



sm 0" 



oo 



and all expansions are uniformly on any compact subset of the set (|1.18p . Let us go back 
to (|16.4|) and (|16.5|) . Since v = uq + h , we notice 

R(s,s) = -iu (8s 2 + 2x) - i(4s 2 + x) 

+ip(16s 2 + Ax) \r' 12 (1)R 22 (1) - ^'22(1)^12(1)] + 0(s~ 1 ) 

and similarly 

R(-s,-s) = -iu (8s 2 + 2x) - i(4s 2 + x) + ip(16s 2 + Ax) 



R' U (-1)R 21 (-1) - R' 2l (-1)R U (-1) + Ois- 1 ) 



Next 



^'l2 


(1) = 


I(Q(1) + (I + 5)Q'(1)) 


R22 


(1) = 


|((/ + S)0(l)) 2a 


R 22 


(1) = 


I(Q(1) + (J + 5)Q'(1)) 


Rl2 


(1) = 


i((/ + fl)Q(l)) 12 



1 12 



'22 



and therefore 



1 r, 



^12(1)^22(1) - ^(1)^12(1) = i [(Q(l) + (1+ B)Q'(1)) 12 ((/ + B)Q(1)) 22 
-(Q(l) + (/ + B)Q'(1)) 22 ((J + B)Q(1)) 12 " . 
Now combine the previously derived information on Q(1),Q' (1) and i? to derive 

ip(l6s 2 + Ax)\r' 12 (1)R 22 {1) -^22(1)^12(1)] = (4s 2 + x)(i + tancr) +a + + 0(s~ 1 ) 
with a function a + = a_|_(,s, 2, 7) such that 

Ja + (s,x, 1 )ds = 0(lns), s^oo. 

Following the same computations for R(—s, —s) 

#(_ S) _ S ) = -iu {8s 2 + 2x) + (4s 2 + x)tan a + + O(s^) . (18.7) 
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where 

and together from Proposition [1] 
d 

— In det(J - Tocsin) = ^o(16s 2 + 4s) - (8s 2 + 2x) tan a - a+ - a_ + 0{s~ l ) . (18.8) 
as v ' 

Opposed to the latter equation we now recall Proposition [5] and evaluate the logarithmic 

^-derivative. For 7 > 1 

Xi = lim (x(X(X)e- i{ ^ i+xX)aa - I) 

A— s>oo V 



-2vsa z + s(a 3 + B) + ^ ^ Q_ (w) (w) (Gq (w) - I) dw 



where the expansion for B has already been computed. From this and residue theorem 
v _ is /sincr — 1 \ i ( v —u 

Ai = — AVSGj, + S0"3 

COS (7 



1 — sincr / 2 U —v 



i(v sincr — u) /sincr — 1 



cos^ a 



2 a \ 1 — sin <T 



)+o( s - 1 ) ; 



hence 



d 2zz sin cr 
— lndet(J - 7K csin ) = 4iv s + v - 2s tan a = 1 5 hO(s~ 1 ). (18.9) 

OX COS^ (7 COS^ C7 

Integrating both identities (|18.8p . (|18.9p and comparing the result, we conclude for 
s — > 00 outside the zero set (|18.6p 

lndet(J — 7-K"csin) = ivo(^s 3 + 4sx^J +ln\cosa(s,x,j)\ + xi\ns 

00 

- J{V- x)u 2 (y, i)dy + X2 + Ois" 1 ) , (18.10) 

X 

with real- valued constants Xh solely depending on 7 and the error term is uniform on any 
compact subset of the set (|1.18p . The given expansion (|18.10p verifies the claim on the 
asymptotic distribution of the zeros of the Fredholm determinant as given in Theorem 
1 

Remark 3. We want to emphasize that our proof of Theorem^ produces in fact an entire 
asymptotic series for det(/ — 7-fC CS m) in case 7 > 1, similarly to what we obtained in 
Theorem [2. However due to the increased amount of computations within the Riemann- 
Hilbert analysis, we chose not to compute the constant Xi(l)- 

REFERENCES 

[1] E. Basor, C. Tracy, Some problems associated with the asymptotics of r-functions, translated into 

Japanese for Swikagaku, 345, (1992) 
[2] E. Basor, H. Widom, Toeplitz and Wiener-Hopf determinants with piecewise continuous symbols, 

Journal of Functional Analysis, 50 (1983), 387 - 413. 
[3] H. Bateman, A. Erdelyi, Higher Transcendental Functions, McGraw-Hill, NY, 1953. 



ASYMPTOTICS OF A CUBIC SINE KERNEL DETERMINANT 



58 



[4] M. Bertola, On the location of poles for the Ablowitz-Segur family of solutions to the second 

Painleve equation, preprint: arXiv: 1203.2988 
[5] B. Bettelheim and P. Wiegmann, Fermi distribution of semiclassical non-equlibirum Fermi states, 

Phys. Rev. B, 085102 (2011) 
[6] A. Budylin, V. Buslaev, Quasiclassical asymptotics of the resolvent of an integral convolution 

operator with a sine kernel on a finite interval, Algebra i Analiz, 7, no. 6 (1995), 79-103 
[7] P. Bleher and A. Its, Semiclassical asymptotics of orthogonal polynomials, Riemann-Hilbert prob- 
lem, and universality in the matrix model, Annals Math.150 (1999), 185-266. 
[8] P. Bleher and A. Its, Double scaling limit in the random matrix model: the Riemann-Hilbert 

approach, Comm. Pure Appl. Math. 56 (2003), 433-516. 
[9] T. Bothner and A. Its, Asymptotics of a Fredholm determinant involving the second Painleve 

transcendent, preprint: larXiv:12 09.5415 (2012) 
[10] T. Bothner and A. Its, The nonlinear steepest descent approach to the singular asymptotics of the 

second Painleve transcendent, Physica D 241 (23-24) (2012) 2204-2225 
[11] T. Claeys and A. Kuijlaars, Universality of the double scaling limit in random matrix models, 

Comm. Pure Appl. Math., 59, no. 11 (2006), 1573-1603. 
[12] P. Deift, Orthogonal Polynomials and Random Matrices: A Riemann- Hilbert Approach, Courant 

Lecture Notes in Mathematics 3, Amer. Math. Soc, Providence, RI, 1999. 
[13] P. Deift, A. Its, I. Krasovsky, On Asymptotics of a Toeplitz determinant with singularities, preprint: 

larXiv: 1206. 12921 (2012) 

[14] P. Deift, A. Its, I. Krasovsky, and X. Zhou, The Widom-Dyson constant for the gap probability in 
rndom matrix theory, J. Comput. Appl. Math. 202 (2007), no. 1, 26-47 

[15] P. Deift, A. Its, and X. Zhou, A Riemann-Hilbert approach to asymptotic problems arising in the 
theory of random matrix models, and also in the theory of integrable statistical mechanics, Ann. 
Math. 146 (1997), 149-235 

[16] P. Deift, T. Kriecherbauer, and K.T-R McLaughlin, New results on the equilibrium measure for 
logarithmic potentials in the presence of an external field, J. Approx. Theory 95 (1998), 388-475 

[17] P. Deift, T. Kriecherbauer, K.T-R McLaughlin, S. Venakides, and X. Zhou, Uniform asymptotics 
for polynomials orthogonal with respect to varying exponential weights and applications to uni- 
versality questions in random matrix theory, Comm. Pure Appl. Math. 52 (1999), 1335-1425 

[18] P.A. Deift and X. Zhou, A steepest descent method for oscillatory Riemann-Hilbert problems. 
Asymptotics for the MKdV equation, Ann. of Math., 137 (1993) 295-368. 

[19] F. Dyson, Fredholm determinants and inverse scattering problems, Commun. Math. Phys. 47 (1976) 
171-183 

[20] T. Ehrhardt, Dyson's constant in the asymptotics of the Fredholm determinant of the sine kernel, 

Comm. Math. Phys. 262 (2006), 317-341 
[21] L.D. Faddeev and L.A. Takthajan, Hamiltonian Methods in the Theory of Solitons, Springer- Verlag, 

Berlin, Heidelberg, 1987 

[22] A. Fokas, A. Its, A. Kapaev, V. Novokshenov, Painleve transcendents. The Riemann-Hilbert ap- 
proach. AMS series: Mathematical surveys and monographs, 128 (2006) 553p. 

[23] H. Flaschka and A.C. Newell, Monodromy- and spectrum-preserving deformations I, Comm. Math. 
Phys. 76 (1980) 65-116. 

[24] A. Its, A. Izergin, V. Korepin, Long-distance asymptotics of temperature correlators of the impen- 
etrable Bose gas, Commun. Math. Phys. 130, 471 (1990) 

[25] A. Its, A. Izergin, V. Korepin, and N. Slavnov, Differential equations for quantum correlation 
functions, Int. J. Mod. Physics B4 (1990), 1003-1037 

[26] A. Its, A. Izergin, V. Korepin, G. Varzugin, Large time and distance asymptotics of the temperature 
field correlator in the impenetrable bose gas, Nucl. Phys. B 348, 757-765 (1991). 

[27] A. Its, A. Izergin, V. Korepin, G. Varzugin, Large time and distance asymptotics of field correlation 
function of impenetrable bosons at finite temperature, Physica D 54, 351-395 (1992). 

[28] A. Its, I. Krasovsky, Hankel determinant and orthogonal polynomials for the Gaussian weight with 
a jump, Contemp. Math. 458 (2008), 215-247 



ASYMPTOTICS OF A CUBIC SINE KERNEL DETERMINANT 



59 



[29] M. Jimbo, T. Miwa, and K. Ueno, Monodromy preserving deformation of linear ordinary differential 
equations with rational coefficients. I. General theory and r-function, Phys. D 2 no. 2 (1981), 306- 
352. 

[30] N. Kitanine, K. Kozlowski, J. Maillet, N. Slavnov, V. Terras, Riemann-Hilbert approach to a 
generalised sine kernel and applications, Comm. Math. Phys. 291 (2009) 691-761 

[31] I. Krasovsky, Gap probability in the spectrum of random matrices and asymptotics of polynomials 
orthogonal on an arc of the unit circle, Int. Math. Res. Not. 2004 (2004), 1249-1272 

[32] M.L. Mehta, Random Matrices, 2nd edition, Academic Press, Boston, 1991. 

[33] B.M. McCoy and Sh. Tang, Connection formulae for Painleve V functions, Physica D 19 (1986), 
42-72 

[34] B.M. McCoy and Sh. Tang, Connection formulae for Painleve V functions: II. The <5-function Bose 

gas problem, Physica D 20 (1986), 187-216 
[35] H. Montgomery, The pair correlation of zeros of the zeta function, Proc. Sympos. Pure Math., vol. 

24, Amer. Math. Soc, Providence, R.I., 1973, pp. 181-193. 
[36] A. Odlyzko, On the distribution of spacings between zeros of the zeta function, Math. Comp., 48 

(1987), pp. 273-308 

[37] L. Pastur, M. Shcherbina, Universality of the local eigenvalue statistics for a class of unitary 
invariant random matrix ensembles, J. Statist. Phys. 86 (1997), 109 - 147. 

[38] H. Segur, M. Ablowitz, Asymptotic solutions of the Korteweg de Vries equation, Stud. Appl. Math. 
571 (1977), 13-44 

[39] B. Simon, Trace ideals and their applications, London Mathematical Society Lecture Note Series, 

35. Cambridge University Press, Cambridge-New York, 1979 
[40] B.I. Suleimanov, On Asymptotics of Regular Solutions for a Special Kind of Painleve V Equation, 

Led. Notes in Math., Springer Verlag, 1191 (1986) 230-255. 

Department of Mathematical Sciences, Indiana University-Purdue University Indianapo- 
lis, 402 N. Blackford St., Indianapolis, IN 46202, U.S.A. 
E-mail address: tbothner@iupui.edu 

Department of Mathematical Sciences, Indiana University-Purdue University Indianapo- 
lis, 402 N. Blackford St., Indianapolis, IN 46202, U.S.A. 
E-mail address: itsa@math.iupui.edu 



